General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



A Final Report 
Grant No. NSG- 1652-2 

TERRESTRIAL GROWTH OF LEAD-TIN-TELLURIDE BY TECHNIQUES 
RELATED TO LOW G GROWTH 


Submitted to: 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665 

Attention: Archibald L. Fripp, Jr. 

FED 

Submitted by: 

W. A. Jesser 
Professor 



Department of Materials Science 
RESEARCH LABORATORIES FOR THE ENGINEERING SCIENCES 
SCHOOL OF ENGINEERING AND APPLIED SCIENCE 
UNIVERSITY OF VIRGINIA 
CHARLOTTESVILLE, VIRGINIA 

(NASA-Ca- 169541) TERRESTRIAL GROWTH Of N83-13230 

L£AD-TIN-TELLURID£ BY TECHNIQUES RELATED TO 
LOW G GROBTd final Report (Vir<jinia Univ.) 

153 p HC A08/HF AOl CSCL Ilf Uncias 

G3/26 02073 


Copy No . 






Report No. UVA/528185/MS82/101 
October 1982 


ORIGINAI. O 
OF POOR QuALiTY 


1. Btport No. 


2. Gov«rnm«nt Acctuion No. 


3. R«cipi«m‘i Ciuiog No. 


4 . Title end Subtitle 

TERRESTRIAL GROWTH OF LEAD-TIN- TELLURIDE BY TECHNIQUES 
RELATED TO LOW G GROWTH 


7. Authof(s) 

W. A. Jesser, Professor 

T. I. Ejim, Graduate Research Assistant 


9. Performing Organization Name and Addresa 

Department of Materials Science 
University of Virginia, Thornton Hall 
Charlottesville, VA 22901 


12. Sponsoring Agency Name and Address 

NASA Langley Research Center 

Hampton, VA 23665 

Attn: Archibald L. Fripp, Jr. , FED 


15. Supplementary Notes 


5. Report Date 

November 1982 


6. Performing Organization Coda 

5-28185 


8. Performing Organization Report No. 

UVA/528185/MS82/101 


10. Work Unit No. 


11. Contractor Grant No. 

NSG-1652-2 


13. Type of Report and Period Covered 
Final Report f 


14. Sponsoring Agency Code 


16. Abstract 

A modified Bridgman-Stockbarger furnace was constructed for a study of the solid- 
ification of silver, germanium and lead-tin-telluride . The melt-solid interface posi- 
tion with respect to the furnace and its temperature profile was determined by mea- 
suring the discontinuity in the slope of temperature as a function of position in the 
melt and in the solid. The results show that the interface position of the semi- 
conductors germanium and lead-tin-telluride was essentially constant with respect to 
the furnace, and hence the growth rate was constant and equal to the sample trans- 
lation rate of 0.046 cra/min and 0.178 cm/min in each case. The metal, silver, on the 
other hand showed a continuous interface migration toward the hot zone of the furnace 
and always exhibited a growth rate which was higher than the ampoule translation rate. 

The K^/Kq ratio of lead-tin-telluride was determined to be 2.33 i 0.06 where K, 

denotes the thermal conductivity of the liquid, solid respectively. The value of K.. 

- 1-1 ^ 
was calculated to be about 0.054 Watt cm K 

The diffusion boundary layer thickness was calculated for lead-tin-telluride to 
be about 0.05 cm using a liquid diffusivity of 7x10 cm^/sec. 


17 Key Words ISuggetTed by Author(s)) 18. Distribution Statement 

Bridgman-Stockbarger furnace, 
lead-tin-telluride, germanium, single 

crystal growth, directional solidif icat ior Unclassified-Unlimited 


19 Security O^uif (of this reporti 

20 Security Claisif (of this page) 

21. No. of Paget 

Unclassified 

Unclassified 

154 



jOS 


Fch sale by Ihe National Technical Infoimation Service. Springfield. Virginia 22161 





















OF POOR 

ACKNOWLEDGEMENTS 

This fiaal report is the Ph.D. dissertation of Dr. Theo Ejim who 
completed his Ph.D. work at the University of Virginia in September 1982 
after spending two years of effort in the NASA Langley Research Center 
in Hampton, Virginia. Extensive and generous use was made of NASA 
equipment which is not available at the University. The guidance and 
advice of NASA staff in the use of equipment and during the pursuit of 
the goals of this research is gratefully acknowledged. 


ii 



TABLE OF CONTENTS 


Page 


Abstract ^ 

Acknowledgements 11 

List of Figures vl 

List of Tables 

List of Symbols 

CHAPTER 


I. lOTRODUCTION 1 

II. BACKGROUND AND PREVIOUS WORK 6 

A. Interface Breakdown 6 

1. Constitutional Supercooling Criterion 6 

2. Stability Criterion 9 

B. Mass Redistribution 10 

1. Liquid Diffusion Controlled 10 

2. Partial and Complete Mixing In the Melt 12 

C. Heat Analysis 13 

1. Effect of Ampoule Motion 13 

2. Tickler Heater Analysis Id 

D. Convection In the Melt 20 

III. APPARATUS AND EXPERIMENTAL DETAILS 25 

A. Apparatus 25 

1. Heater Blocks 25 

2. Controllers and Isothermal Liners 27 


iii 



CHAPTER 


Page 

3. Tickler Heater 28 

4. Ampoule Design and Vacuum System 22 

B. Experimental Technique 34 

1. Bridgman Technique 34 

2. X-ray Analysis 35 

C. Experimental Procedure 37 

1. Furnace Characterlzatloii ' 37 

2. Melt-solid Tempera <'ure Measurement 39 

3. Quantitative X-ray Analysis 44 

D. Sample Preparation '44 

1. C''mpoundlng Pb Sn Te 44 

1"X X 

2. Cleaning and Etching of Samples 45 

3. Samples for X-ray Analysis 47 

IV. EXPERIMENTAL DATA AND RESULTS '....... 49 

A. System Geometry 49 

B. Furnace Characterization 49 

C. Melt-solid Temperature Measurement 63 

1. Interface Position Detemination 63 

2. Materials Constant 67 

3. Tickler Heater Effects 70 

>4. Translation Rate Effects 75 

5. Grovth Rate Variations 78 

B,. Composi tlon Analysis 

V. ERROR ANALYSIS 92 

A. Melt-solid Thermal Coiductivity Ratios 92 

B. Position Error 103 

C. Concentration Profile Error 105 

iv 



CHAPTER 


Page 

VI. DISCUSSION 106 

A. Furnace Characterization 106 

1 . 'i.ntroductlon 106 

2. Effect of Plugging 107 

. 3. Effect of Loading 108 

4. Tickler Heater Effect 110 

B. Melt-solid Interface Positions llO 

1. Introduction 110 

2. Tickler Heater Effect on Interface Position 112 

3. Effect of Translation Rate 113 

4. Materials Constant 114 

C. Growth Rate 118 

D. Solute Redistribution 120 

VII. SUMMARY AND CONCLUSIONS ..... 125 

A. Growth Rates 125 

B. Materials Constants 125 

C. Plugging and Tickler Heater Effects 127 

D. Load Material and Translation Rate Effects 128 

APPENDIX I - Data Reduction Program 130 

APPENDIX II - Examples of Interface Position and Freezing Poin . 

Determination Curves 133 


BIBLIOGRAPHY 


142 


LIST OF FIGURES 


Figure Page 

1. Condition for Constitutional Superheating 23 

2. A Two Zone Br:^dguan Furnace 26 

3a. Photograph of the Constructed Furnace 29 

3b. Schematic of the Final Form of the Furnace 30 

4. The Constructed Tickler Heater 31 

5. Ampoule Design Used for Growth 33 

o . The Vacuum System 40 

7. Directionally Solidified Rods of Pb Sn Te and Ag 43 

JL*X X 

8a. PbTe Phase Diagram 46 

8b. Pb. Sn Te Pseudobinary Phase Diagram 46 

X"X X 

8c. Ternary Phase Diagram of PbSnTe 46 

9. System Geometry and Definition of Coordinates 50 

10. Unplugged Empty Furnace Profile 52 

11. Plugged Empty Furnace Profile 53 

12. Isotherms of a Plugged Empty Furnace 55 

13. Plugged Empty Furnace Profile With 1.43 cm Long 

Tickler Heater 56 

14. Temperature Gradient vs. Distance Along the Furnace 57 

15. Plugged Empty Furnace Profile With a 2.3 cm Long 

Tickler Heater . 58 

16. A Silver Loaded Furnace Profile 60 

17a. Empty Amooule Loaded Furnace Profile 61 

17b. Isotheims o2 the Loaded Furnace 62 

18. Temperature vs. Distance From Melt Into the Sc’ id 64 

vi 



Figure Page 

19. Temperature Gradient vs. Distance From Melt 

into the Solid 65 

20. £:cpanded Form of Figure 18 and 19 66 

21. Effect of Ticklti Heater Length on Interface Positions 72 

22. Effect of Tickler Heater Power Setting on Interface 

Positions 74 

23. Translation Rate Effect on Interface Position 77 

24. Shift ir. Temperature Profile With Freezing 79 

25. Growth Rate Changes 81 

26. Fraction Solidified vs. Translated Length 82 

27. Concentration Profile 83 

28. PbSnTe X-ray Spectrum 86 

29. X-Ray Spectrum From all Etched Grain Boundaries 91 

30. Parallel Resistor Analog> 94 

31. Curvature Effect 101 

32. Possible Heat and Mass Flow Pattern for the Indicated 

Interface Shapes 101 


vii 



ORIGI!«/vL i 
OF POOR 


LIST OF TABLES 


Table 


1.. Effect of Lead Material on Temperature Gradient 69 

2. Comparison of Measured Melt-solid Thermal Conductivity 

Ratios With Selected Literature Values 71 

3. Effect of Tickler Heater Length on Temperature Gradient 

Ahead of the Interface 73 

4. Measured Compositions Showing a One-to-one Ratio of 

Metal and Nonmetal 88 

5. X-ray Diffraction Data Measured 89 

6. X-ray Diffraction Data From ASTM Cards 89 

7. Concentration Versus Temperature Corrected for 

Interface Curvature 123 


I 

i 

1 

I 


'% 

4e 

4 

1 


viii 


LIST OF SYMBOLS 


Cross sectional area 

Concentration of element "A" in the sample 
Solute concentration in the melt 
Heat capacity of constant pressure 
Solute concentration in the solid 
Solute diffusion coefficient in the melt 
Solute diffusion coefficient in the solid 

Position of thermocouple junction with resy>ect to the bottom 
of the ampoule 
Volume fraction 

Concentration gradient in the solid 

Temperature gradient in the melt 

Temperature gradient in the solid 

fraction solidified 

Latent heat of fusion 

Heat transfer coefficient in the melt 

Heat transfer coefficient in the solid 

Interface position with respect to FRP 

Normalized interface position, (Z - i)/L^ 

X-ray intensity due to pure element "A" 

X-ray intensity due to element "A" in the sample 

Thermal conductivity in the melt 

Thermal conductivity in the solid 

Solute segretation coefficient, C_/C 

^ L 

Length solidified at a given instance, Z - i 



Total length of grown rod (boule) 

Slope of the liquidus curve versus composition 
m_ Slope of the solidus curve versus composition 

Pe Peclet number 

q Heat flow rate 

R Growth rate 

5 Radius of the sample cross section 

t Time 

T Equilibrium liquidus temperature 

Li 

TC Thermocouple 

TH Tickler heater 

T^ Cold zone temperature 

T Hot zone temperature 

V Translation rate 

X' Distance from the interface into the melt 

Z Translated length measured from the FRP 

Bottom of ampoule with respect to FRP 
a Length of tickler heater 

6 Blot number 

p Density 

5 Width of diffusion boundary layer 


oRiGirj'L F is 

OF POOR QUALITY 


X 



CHAPTER I 


INTRODUCTION 

Directional Solidification (DS) and normal freezing carried out 

in a controlled environment are particularly important in the field 

of single crystal growth and, therefore, to the electronics Industry. 

The Bridgman-Stockbarger^^^ technique is one of the common methods 

of growing bulk single crystals. This technique draws heavily from 

DS studies. The one common aim of crystal growers is to produce a 

homogeneous sinple crystal of contj oiled composition free of line, 

point and volume defects, which are detrimental to device applications 

of the crystals. Some important controllable factors which determine 

the extent to which the grower's aim is met include the rate of 

growth and the temperature gradients in the phases. These factors 

( 2 ) 

Influence the heat flow pattern in a given material. 

In the Bridgman technique used in the present work, much like DS 
experiments, heat is Introduced into the material from a heater, through 
an adiabatic zone and out of the cc oier in a vertical mode. The cooling 
device in a typical system can be coolant such as water or liquid metal 
for direct immersion, or in the present case it can be a heater block 
with a lower temperature setting than the hot zone. The choice of the 
cooling method may depend on the knowledge of the temperature gradient 
necessary for the avoidance of constitutional supercooling, the magnitude 
of the thermal gradient ana cooling rates that will not introduce defects 
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in the grown crystal by the subsequent thermal strain and point 
defect asgregation. 

The macroscopic melt-solid interface shape may greatly influence 
the generation and '"■opagation of defects during directional solidifi- 
cation. It is generally accepted that a convex shape (looking from 
the melt side) favors grain selection. ' The well known constitu- 
tional supercooling concept indicates that slow rates of growth 
and c high temperature -gradient in the liquid phase are necessary 
cbndltions for homogeneous ;wth of alloy materials. These concerns 
are particularly important in the field of semiconductor device 
applications. Developments in this field call for very high quality 
single crystals in order to achieve uniform electrical response within 
the crystal. This means, among other things, the homogeneous distri- 
bution of dopants and solutes and a very low dislocation density. In 
an alloy semiconductor like Pb. Sn Te, homogeneity means the absence 
of convection in the melt end the avoidance of constitutional super- 
cc'jling with its resultant dendritic growth. The absence of convection 
can only be achieved in a system that exhibits both thermal and solutal 
stability in the melt at all ffmes during the growth. If for any con- 
figuration, the lighter material (lighter either by differences in 
temperature or solute concentration) is at the bottom of the system 
then the density gradient can couple with gravity to start convection in 
the melt. Looking at the phase diagram of binary alloy system and 
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employing the lever rule, it is easily seen that the composition of 
the first-to-freeze solid may be higher or lower than the bulk melt 
composition depending on the appropriate segregation coefficient. 

The segregation coefficient, here, means the ratio of the solute 
coacentration in the solid phase and the solute concentration in 
the melt from which the solid froze. This rejection or depletion 
of solute creates a diffusion boundary layer in the melt Imme*. lately 
ahead of the interface. Within this layer mass transfer should be 
by concentration gradient-driven diffusion alone resulting in homo- 
geneous steady-state growth. However, this is not always true and 
much of the quality of the grown crystal depends on the behavior of 
this layer. If any region of this layer sees a temperature that is 
lower than the freezing point, that region will break down leading to 
dendritic or cellular growth. Also strong gravity-driven convection 
can sweep away this boundary layer into the bulk melt leading to 
nonsteady state and nonhomogeneous growth. 

The microgravity environment of space is being proposed as an 
environment in which the problem associated with gravity-driven con- 
vection can be lessened if not entirely eliminated. The present 
research is part of a general effort to develop experiments for materials 
processing in space aboard the space shuttle. It is, tn refore, a part 
of an Intensive attempt to understand all factors that can influence 
crystal quality. Until recently the ampoule translation rate during 
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directional solidification had been regarded as the crystal growth 
rate. The two rates are now known to be able to differ by as much 
as lOOZ.^^^ Experimental furnace characterization has been minimal. 
While there has been much theoretical work in this area, the mathe- 
matical difficulties Involved have led to the introduction of 


assumptions that may not fit experimental conditions. Much experi- 
mentation is still needed in order to understand the thermal environ- 
mental changes that take place in a Bridgman-^ tockbarger system during 
crystal growth. 

Three different materials were used in this investigation. These 
include two semiconductors, germanium and Pbj^^^Sn^Te and one metal, 
silver. The alloy semiconductor Pbj^ ^Sn^Te is a technologically very 
important material being proposed as one of the materials to be pro- 
cessed in the space laboratory experiments. It has an energy gap that 
is dependent on the composition and can be as low as zero.^^^’^^^ 

With the low energy gap the advantage that its energy gap can be 
controlled by varying lomposition, the semiconductor may be used in 
photovoltaic detectors and laser over a range of wavelengths. 

It has been mentioned earlier in this chapter that the shape of 

( 3 ) 

the melt-solid interface influences grain selection. Theoretical 

( ■?) ( l) 

work also indicates that the position of the melt-solid inter- 

face with respect to the furnace determines the shape of the interface. 


\ 


i 

» 


The present work concentrates on understanding how the melt-solid 
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interface position changes in relation to a furnace reference point 
and with respect to a number of furnace parameters and ampoule trans- 
lation rates. With the knowledge of the interface movements the 
actual growth rate of the material will be determined as opposed to 
the translation rate. A careful characterization of the empty furnace 
will be carried out in order to understand possible shifts in temperature 
profile on loading. Finally, the concentration profile of a grown 
alloy will be correlated with the freezing temperature determined at 
selected points along the crystal. 
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CHAPTER II 

BACKGROUND AND PREVIOUS WORK 

A. INTERFACE BREAKDOWN 

1. Constitutional Supercooling Criterion 

One of the many published works dealing with the prrblems of 

.. - : ( 4 ) 

solidification is the classic paper by Rutter and Chalmers in 

which they treated the problem of Constitutional supercooling (CS) 
during the directional solidification of single crystals of tin 
doped with lead. For a given set of growth conditions like dopant 
concentration and temperature gradient in the melt, a critical growth 
rate was reached at which interface breakdown occurred leading to a 
hexagonal array of protrusions. They explained this "cellular" 
structure with the theory of constitutional supercooling. This 
supercooling ahead of the interface is a result of either solute 
buildup or depletion in the melt immediately ahead of the liquid- 
solid interface. Factors that affected the interface breakdown 
included the rate of translation of the ampoule or withdrawal of the 
furnace, the amount of impurity and the temperature gradient in the 
melt. To develop quantitatively the C.S. criterion, the heat and 
mass flow are considered only at the liquid-solid interface. The 
gradient of the solute in the liquid at the interface for a steady- 
state condition is given as 
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( 1 ) 


where « solute concentration in the p.elt at the interface 

R » growth rate 

■ solute diffusion coefficient in the melt 
k * solute segregation coefficient, C„/C* 

D Xrf 

Cg « solute concentration in the solid 
X' ■ distance from the interface into the melt 
Assuming equilibrium at a flat interface the derivative of the 
equilibrium liquidus temperature T with respect to distance from 

Li 

the interface into the melt is given as 



where ™ “ the slope of the liquidus curve versus composition. Con- 

stitutional supercooling is absent when the actual temperature gradient 
in the melt is equal to or greater than the gradient of the liquidus 
temperature curve. 

Assuming that k is independent of composition and substituting 

/Q\ 

Eq. (1) into Eq. (2) gives the C.S. criterion^ for interface breakdown. 
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^ K|C*3(l-k) 

R ^ Dj^k 


( 3 ) 


where ■ actual temperature gradient in the melt 

L 

C* ■ solute concentration in the solid at the interface 
Hurle^^^ has treated the case where the melt is stirred. The tem- 
perature gradient in the melt at the interface is given as 





Im^lR 




(4a) 


where C* ■ the steady-state melt concentration at the interface 
C«j ■ solute concentration in the solid 
He obtained the C.S. criterion for interface breakdown with a stirred 
melt as 

G |m |C (1-k) 

_L < L 2 ••• (4b) 

R — k+(l-k) exp(-R<5/Dj^) 

where C “ average composition of the bulk melt 
o 

5 - width of a diffusion layer 

In deriving Eq. (4b) he used the effective distribution coefficient 
proposed by Burton et for a system with some diffusion layer 

of width 6, within which mass transport is by diffusion only. This 
will be discussed later in this chapter. 
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2. Stability Criterion 

The C.S. criterion of Eq. (3) is basod on static thermodynamic 
considerations. A supercooled region is one where the solid phase of 
the material is thermodynamically more stable than the melt and there- 
fore crystals at the interface with the right orientation will rapidly 
grow out to this region, leading to an interface breakdown. 

A more fundamental theory that describes interface breakdown was 
developed by Mullins and Sekerka. The Interface stability theory 
differs from the C.S. principles in two main respects. The theory is 
based upon the dynamics of the whole system rather than upon the 
thermodynamics of the melt ahead of the interface. Such things as 
heat flow in the solid and latent heat of fusion are not accounted 
for in the C.S. principle. The C.S. principle deals with the question 
of which state, liquid or solid is thermodynamically favorable in the 
region of melt ahead of the Interface, but the real question is what 
state is dynamically achievable. Second, the stability theory of 
Mullins and Sekerka^^^^ describes the time evolution of a perturbed 
interface and the resulting concentration and temperature fields. 

With the assumption that the bulk of the phase was isotropic and that 
a local equilibrium existed at the interface, they solved the Laplace 
equations of both the thermal and mass fields taking into account the 
solid-liquid interfacial energy, latent heat of fusion and the thermal 
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conductivities of both the liquid and solid phases. They conclude 
that high surface tension and a high temperature gradient favor inter- 
face stability whereas the solute buildup or depletion, depending on 
the value of the segregation coefficient, is destabilizing. 

Ignoring the "capillary*' term (surface tension) which is always 
stabilizing, the interlace stability criterion can be simplified to the 
form below. 






^ \ 
H 


R < 0 


(5) 


where K_ . ■ thermal conductivity in solid, liquid 

w } JL 

H « latent heat of fusion 
G_ ■ temperature gradient in the liquid 
G_ ■ solute concentration gradient in the liquid 
The latent heat term makes the interface more stable ’■han would be 
predicted by the C.S. criterion. Again, if >> Kg the conductivity 
term makes the interface more stable than predicted by the C.S. criterion. 
It is therefore possible to have constitutional supercooling without 
an instability and instability without constitutional supercooling. 

B. MASS REDISTRIBUTION 

1. Liquid Diffusion Controlled 

A quantitative treatment of the solute redistribution and the 

/3) (12) 

C. S. criterion was developed by Tiller et al. ' and Smith et al. 
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The segregation of Impurities In a dilute binary alloy system Is 
descrlbedby the segregation coefficient. This Is the ratio of the 
solute concentration In the solid to the concentration In the melt 
from which the solid forms. They assumed a system with a segregation 
coefficient of less than unity and negligible convective mixing In 
the liquid. With the above assumptions, they solved the following 
diffusion equation 




3t 


( 6 ) 


where t ■ time. 

Solving the time Independent form of the equation gives the solute 
concentration of the frozen portion and the liquid ahead of the Inter- 
face respectively as 


Cg - Cq {(1-k) [l-exp(-k 1^)] + k} (7a) 

L 

=L- S 


where Is the starting bulk melt concentration. 
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2. Partial and Complete Mixing In the Melt 
In most cases of crystal growth from the melt, convection Is 
present and solute redistribution In these cases has been treated by 
Burton et al.^^^^ and by Wagner. In a solidification process 
from the melt of an alloy, there Is always a solute buildup or depletion 
In the melt Immediately ahead of the solid-liquid interface. In their 
treatment. It Is assumed that t^iere Is a diffusion layer of width 6, 
which Is Independent of the convection In the bulk melt. Outside this 
layer, the melt composition is maintained uniform by convection and 
Inside it, mass transport is by diffusion only. Again the same time 
Independent diffusion equation is solved but with modified boundary 
conditions: 


(x*"0) at the interface 
C ■ C_ (x’^6) in the bulk melt 

Li U 

The solution to the diffusion equation is 



(8a) 

(8b) 


(9) 


where C *, C * are liquid and solid solute concentrations at the 

la 3 

interface 

6 ■ width of the diffusion boundary layer 
- average concentration in the bulk melt 
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In the same work, they proposed the value of an effective 
segregation coefficient, as 


eff 


-R6/D, 


( 10 ) 


k +(l-k ) e 


where k • equilibrium segregation coefficient, C^/C^* 


Kti ■ 


( 14 ) 


As indicated by Flemings, ■ besides the treatments of mass 
redistribution by Tiller et al.^®^ and by Burton et al.^^^^ discussed 
above, there is the earlier treatment by a number of invest!- 
gators^^^^ ’ in which complete mixing in the liquid is assumed. 
Consider a column of melt of an alloy material being lowered into the 
cold zone of a furnace in a normal freezing process. The first to 
freeze is kC^. As cooling and solidification progresses, there is 
a solute buildup in the melt near the interface. At temperature T* 
(temperature at the interface), a solid of composition Cg* freezes 
out of a liquid of composition Equating the solute rejected 

when a small amount of solid forms with the resulting solute increase 
in the melt will give 


* Cg*) dg^ . (1-g^) dCj^ 


( 11 ) 


Substituting the equilibrium segregation coefficient and integrating 
from C-,* ” kC at g ■ 0 yields the composition of the solid at 

bo S' 
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the liquid-solid interface as 

S * ■ '“=0 


( 12 ) 


where C ■ average composition in the melt, 
o 


C. HEAT ANALYSIS 


1. Effect of Ampoule Motion and Insulation Thickness 

(3) 

A theoretical heat analysis work by Chang et al. has con- 
sidered the effects of ampoule movement (Feclet number) on the shape 
and position of the isotherms, in a Bridgman- type furnace. The work 
assumes an infinitely long ampoule, isotropic properties (the section 
of the boule in the hot zone and the section in the cold zone have the 
same physical constants) , negligible latent heat of fusion and a step 
function in Che temperature profile of Che furnace. Using these 
assumptions, they solved a two-dimensional steady-state differencial 
equation for the temperature field as given by 


_K_ A 1, ^ 3^ V II 

^Sar^ az2 


0 


(13) 


where K ■ thermal conductivity of the boule 
p ■ density of the boule 
Cp ■ heat capacity of the boule 
V ■ boule translation rate 
r ■ radial position 


z ■ axial position 



The boundary conditions were as follows: 


at 2 -► a», T -► T (cold end temperature) (lAa) 

c 

at 2 T ■> (hot end temperature) (l^b) 

at r ■ 0 (centerline, 3T/3r ■ 0 (symmetry) (l^c) 

at r ■ S (surface of sample) for 2 > 0, -K 3T/3r ■ h(T-T^) (14d) 

at r - S for 2 < 0, -K 3T/3r - h(T-Tj^^ (lAe) 

at 2 - 0, - T_ and (3T/32)_^ - (3T/32)_ (14f) 


where h is the heat transfer coefficient at the sample surface. 

( 18) 

The soluticii to these equations was obtained by Chang' as 


T-T 

c 

T. -T 
h c 


a r 


■ Vo If > 


(15«) 


for 2 > 0 (cold 2 one) 
and 




1 


+ 



c J 
n o 



(15b) 


for 2 < 0 (hot 2 one) 

where J (a r/S) is the Bessel function of the first kind of order 
0 n 

rero of the argument 
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Pe - V CpS/K (16) 

is the Peclet number (ratio of heat carried by motion of sample to the 

heat by conduction) , a are the eigenvalues determined by a J, (a ) - 

n ” n 1 n 

6j (a ) "0 where J, is the Bessel function of the first Kind of 
on 1 

first order. 
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hS 

K 


(17) 


is the Biot number (ratio of heat loss from the sampK surface to the 
heat conduction) . 


C 

m 


e[pe+(Pe^ >♦- 4 a ^) ] 


(6^+a ^) J (a ) (Pe^ + 4 a ^) 
non n 


C 

n 


6[Pe-(Pe^ + 4 a^^) J 


(6^+a ^) J (a ) (Pe^ + 4 a ^) 
non n 


The work showed that the sensitivity of the interface position (a 
chosen isotherm) to relative heater and coolet temperatures is greater 
for small Biot numbers; that means small ampoule dlamoters, ineffective 
heat transfer from the ampoule surface and high thermal conductivity. 
They showed that the interface becomes increasingly concave (looking 
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from the melt) as the ampoule translation rate Inrreases for constant 
thermal conditions. Of course this statement assumes that the inter- 
face does not break down due to the increased growth rate. Three 
factors combine to cause this increasing concavity of the interface. 

These are the decrease in the interface temperature due to segregation 
(alloy system), the sensible h'iat carried by the ampoule, and the latent 
heat liberation at the Interface during freezing. They also pointed out 
that stirrirg of the melt has roughly the same Influence on Che interface 
position and shape as does a greater thermal conductivity in the melt 
than in Che solid. Both move Che interface toward Che cooler and make 

the Interface less convex. 

( 2 ) 

Fu and Wllco^ have theoretically investigated the influence of 
insulation on the stability of the interface shape and position. They 
used a finite difference technique to solve the steady-state two 
dimensional differential equation for heat conduction in a cylindrical 
rod moving with a known axial velocity. Inserting a layer of insulation 
between the heater and cooler strongly decreased the curvature of the 
Isotherms in this region. Thus the axial temperature gradient at the 
planar interface is expected to be uniform over a larger region with 
increasing thickness of the insulation zone. They showed that for a 
stationary ampoule, the planar isotherm lies in the lower part of the 
insulation region when the Biot number (hS/K) in the hot zone is larger 
Chan that in the cold zone. For equal thermal conductivity in both the 
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hot and cold zones, this result means that the rate of heat transfer 
from the furnace wall to the surface of the sample in the hot zone is 
higher than the rate of heat loss. Finally the work showed a shift in 
isotherms toward the cold zone with Increasing Peclet number (increasing 
ampoule velocity) . This simply means that more hert is being pushed 
down toward the cold zone. 

2. Tickler Heater Analysis 

A one-dimensional analysis of a booster heater, otherwise known 
as a "tickler” heater, inserted in the adiabatic zone, has been performed 
by Naumann.^^^^ He assumed that for a poor conductor, characterized by 
high Biot number, the tickler heater would probably introduce a high 
radial gradient which would distort the isotherms near the solid-liquid 
interface; and no appreciable increase in axial gradient without over- 
heating the material near the sample surface. So the analysis was for 

-2 

good conductors with a Biot number of about 10 in which the radial 
gradients may be neglected. With the further assumption of steady 
state, the problem reduces to the classical one dimensional steady- 
state heat flow in a rod. In dimensionless form it is as follows 
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Pe • Peclet number 
e radius of Che cylindrical rod 
Z * axial position 

The analytical solution to the equation may be described as follows: 
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(19) 


where A, B, C» D are fourler coefficients 

w - 

8 u "Biot number in hoc, cold zones and the furnace is 
H y C 

described as a hot zone, from - * < Z ' ” ^3 

tickler heater, from - Z^ < Z < - 

adiabatic zone, B ■ 0 from - Z^ < Z ^ Zj^ 

cold zone, T from Z- < z < ■» 
c 1 

The function of the tickler heater is to compensate for the roll 
off in sample temperature due to the conduction of heat thereby increas- 
ing the temperature gradient in the melt. However it is necessary to 
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restrict the temperature of the tickler heater region such that the 

maximum temperature attainable In any part of the furnace Is T^. This 

Is required In order to avoid a "hump" In the temperature profile since 

such a condition may enhance convection In the melt which in turn intro- 

(24) 

duces Inhomogenelty Into the crystal. Naumann calculated that the 
maximum tickler heater temperature for his furnace configuration which 
Is given as follows: 


(T^-Tg ) 

max 


(1+BZj^) e^®^ + (1-BZ^) e"^®^ - 2 


( 20 ) 


where L « half the length of the tickler heater. 

D. CONVECTION IN THE MELT 

(19) 

Kim et al. have Investigated the effects of destabilizing 
vertical thermal gradients on crystal growth and segregation in an 
inverted (cold zone on top) Bridgman configuration using a Ga-doped 
germanium. It was found that as solidification progressed the melt 
exhibited successively turbulent convection, oscillatory thermal 
instabilities and finally thermal stability. Growth under turbulent 
convection was shown to be characterized by pronounced and rapid 
fluctuations of growth rate. 

The results indicated that the convective melt flow behavior depends 
strongly on the melt height (aspect ratio) and is in principle analogous 
to the result obtained with an InSb system of similar cor .iguration. 
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In these works, turbulent convection was Indicated between the aspect 
ratio of 2.4 and 1.8 In the germanium system and between 3 and 2 In the 
InSb system. As the melt height decreased, the turbulent convection 
gave way to oscillatory thermal Instability and at an aspect ratio of 
about 0.9, thermal stability was Indicated by the results of the 
germanium system. 

( 21 ) 

Recent work by A. I . Witt e^ al. has treated segregation In a 
Bridgman- type furnace configuration <^or crystals grown In the low 
gravity of space. It was established that the segregation behavior Is 
diffusion controlled without Interference from residual gravity, and 
surface tension gradients. This experiment was conducted during the 
Apollo-Soyuz mlsslcii as a part of NASA's program to determine the 
potential of zero-gravity condltlcuf? for basic research and materials 
processing applications. Ga-doped getrmanlum crystals grown on Earth 
were partially remelted and regrown In space while current pulses were 
applied to achieve Interfaclal demarcation. A number of ground-bared 
control experiments were carried out In an Identical furnace to permit 
comparison of growth and segregation behavior in space and on Earth. It 
was shown that crystals from the ground-based experiments exhibited com- 
plete loss of single crystallinity aftei' a growth length of about 2 cm. 
In contrast, single crystallinity was preserved In the space-grown 
crystal up to the last stages of growth in spite of internal matrix 
breakdown due to constitutional supercooling effects which were expected 
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under the prevailing growth conditions. The microscopic growth rate 
on Earth was found to be. within the experimental error, the same as 
that observed in space. Thus, while unavoidable convection in the 
melt on Earth interfers with the segregation process, it had no signifi- 
cant effect on the microscopic crystal growth rate of a doped element. 
Given the remelting process needed in an experiment like the one 

mentioned above, arid some flrst-to-freeze problems, there Is an increased 

( 22 ) 

effort in the investigation of the melt-bacx process. Yue and Clark 
in their work with a magnesium alloy of copper reported a disproportion- 
ately high solute (copper) concentration in the f irst-to-freeze portion 
of the rod. This la probably due to the nucleation problems that arise 
because of the undercooling required to overcome the nucleation energy 
barrier. It is hoped that with the Introduction of a solid-liquid 
Interface this initial transient may be avoided. It is of interest 
therefore to study the meltback process. Of prime importance is the 

stability of the phase boundary during melt-back. The condition of the 

(23) 

interface may influence the crystal structure on regrowth. Woodruff 
has treated the melting process of an alloy following the same method 
used by Mullins and Sekerka,^^^^ but with the additional consideration 
of solute diffusion in the growing phase. It is shown as in figure 1 
that a region of cDnstitutional superheating will occur ahead of the 
melting interface if 




solute segregation coefficient 

cone, gradient in the solid 

temperature gradient in the solid 

slope of solldius curve versus composition 
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CHAPTER III 

APPARATUS AND EXPERIMENTAL DETAILS 


A. APPARATUS 

1. Heater Blocks 

The furnace that was constructed and characterized in this investi- 
gation has a Bridgraan-type configuration consisting of a hot zone, an 
insulation (adiabatic) zone, and a cold zone. The furnace was required 
to withstand temperatures up to llOO'^ C and exhibit a temperature 
gradient of at least 100° C/cm or more at mperatures of not less than 
900° C. The heater blocks were obtained commercially and according to 

manufacturer's specifications, they could withstand temperatures up to 
0 

1400 C. The heater blocks consist of spirally-wound Kanthal resistance 
wires embedded in ceramic cement. These come in halves, two of which 
will match to form a 4- to 12-inch heater with a cylindrical opening 
between them. The first furnace that was constructed and used for some 
of this investigation is shown in figure 2. It consists of the three 
different regions mentioned above. The adiabatic or insulation zone is 
a region between the hot and cold zones within which there are no heater 
wires and therefore no heat generation. This region is well insulated 
with flberfax material (alumina fibers in a binder), to avoid heat loss. 
Both the hot and cold zones are lined with quartz tubes. This is intended 
to protect the heater coils and the isothermal liners from contamination 
in case of ampoule breakage during an experiment. 
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2. Controllers and Isothermal Liners 

Temperature stability is very important in any crystal growth 
system. It is particularly Important in solidification from the melt 
of alloys since temperature fluctuations can introduce compositional 
variations and may lead to cellular growth. The three furnace thermo- 
couples shown in figure 2 are connected to a Eurotherm proportional 
plus integral (reset) temperature controllers. The controller eliminates 
"droop" by moving the proportional band up or down the temperature scale 
until the point at which control actually occurs coincides with the set 
point. "Droop" here means the difference between the set point and the 
actual point of control. A proportional band Is a region around the set 
point within which power cycles on and off automatically without a 
temperature change. With the Eurotherm controllers, temperature oscillv. 
can be eliminated using the appropriate rate and integral control settings. 
All the control thermocouples used in this work were platlnum/rhodlum- 
platinum (r-type) . 

A number of modifications were done on the furnace to get the final 
form used for most of the experiments. This final form is shown in 
figure 3. The furnace has an isothermal linear in both the hot and cold 
zones. It also contains a booster heater which will be referred to as a 
tickler heater from here on. 

The isothermal liners are basically cylindrical heat pipes. The 

( 21) 

inner surfaces of a heat pipe are lined up with a porous capillary wick. 
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This wick is saturated with the liquid phase of a "working" fluid and 
the remaining volume of the tube contains the vapor phase. Heat applied 
at the evaporator (liquid phase region) by an external source is absorbed 
through the vaporization of the working fluid resulting to a pressure 
differential along the heat pipe column. This pressure differential 
drives vapor from the evaporator (heated region) to the condenser 
(colder region) where it condenses thereby releasing the latent heat of 
vaporization. Through this process, heat is transferred from the 
evaporator region to the condenser region. The main function of the 
isothermal liner is to distribute heat more uniformlly all through a 
given region by the wlcking action. Another function of the isothermal 
liner is to add temperature stability to the furnace through its thermal 
mass. 

3. Tickler Heater 

The final form of the furnace constructed and used in this investi- 
gation is shown in the photograph of figure 3a. A schematic of this 
furnace showing all its Important features is figure 3b. One of the 
important features of this furnace is the presence of a tickler heater 
which as a unit is shown in figure 4. A Kenthal A-1 wire which is 
0.081 cm in diameter is wound on a quartz tube of 1.7 cm ID and 2.1 cm 
OD. The wire ■'.£ wound around the quartz tube In such a way that the 
possible magnetic field created by current in one set of coils is can- 
celled by an opposite field created by current in another set of coils. 


29 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPK 



Figure 3a: The Bridgman-Stockbarqer furnace constructed and used in 

this investigation. 
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In effect a non-inductive coil is created by this type of winding. Three 
tickler heaters were made with 16, 14, and 10 turns of wires, respectively. 

i 

These gave heater coil lengths of 2.54 cm, 2.3 cm, and 1.43 cm, respec- 
tively. These units are then introduced into the adiabatic zone as 
the need arose. The introduction of a tickler heater in each case 
reduced the adiabatic zone length (insulation thickness) from 4.45 cm 
without a tickler heater to (4.45 - a) cm were a is any of the 
tickler heater lengths given above. Each turn of the resistance wires 
is insulated from the others to prevent short circuits. This is accomp- 
lished by covering them with a ceramic material. 

4. Ampoule Design and Vacuum System 

Using the furnace described above, initial attempts were made to 
measure the temperature profile in the melt of Pb,_ Sn Te during normal 
freezing, by dipping a thermocouple junction directly into the melt. 

Several coatings were used to try to protect the thermocouple junction 
from being attacked by the melt. However both bare and coated thermo- 
couple junctionj dissolved in the melt during the experiment. This 
problem coupled with the need to measure the temperature profile both 
in the growing solid phase and the melt called :or a special ampoule 
design. An ampoule was designed with a 1-mm walled capillary tube 
attached to protect the thermocouples. This is shown in figure 5. 

In order to grow materials, it was necessary to clean both the 
ampoule and the charge material and finally to pump on the loaded ampoule 
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under a high vacuum. This of course will help evaporate some contami- 
nants that might be in the system. The vacuum unit used in this 
investigation has a cryct;enic pump which acts as the roughing pump and an 
ion pump which has the ulti ate low pressure capability of 10 torr. 
Attached to the vacuum system is a mass spectrometer unit for the 
analysis of the gases in the system as the loaded ampoule out gases. 

Given the high purity ' ^qulrerients in the single crystal growth of 
electronic materialc>, this vacuum system has ax', obvious advantage over 
a mechanical oil s al i pump which might introduce oil droplets or 
vapor into the sample. 

B. EXPERIMENTAL TECHNIQUE 

] . Bridgman Technique 

Generally crystal growth methods can be divided into three cate- 
gories, namely solid-solid, liquid-solid, and vapor-solid transfor- 
mations. In this Investigation the method employed is the liquid-solid 
transformation. A basic consideration in crystal growth by this method 
is the containment of the molten phase. The technique with which con- 
tainment is achieved can have an appreciable effect on the growth 
conditions and the structural and chemical properties of the resulting 
crystal. Specifically, such things as strain, contamination, spurious 
nucleation and chemical reaction may be a direct result of method of 
containment. The melt-solid transformation as a method of crystal 
growth can be divided into three techniques with reference to the type 
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of contaliunent. These sxr.: (a) the melt and solid phases are in 

contact with a container, (b) only the melt Is In contact with a 
container, and (c) no phase Is In a container. These techniques are 
known as Brldgman-Stockbarger, Czochralskl, and floating-zone melting, 
respectively. 

The Brldgman-Stockbarger technique was used In this Investigation. 

A completely molten phase of the sample and the container (ampoule) are 

pulled from the hot zone of a two zone furnace towards the cold zone 

until all the melt column solidifies. There are a number of factors 

worth mentioning here, that pertain to the growth of crystals by this 

technique. An attempt Is made to achieve a planar melt-solid Interface 

( 22 ) ( 2 “)) 

In order to reduce dislocation density.' ^ An attempt Is also 

made to eliminate or minimize nonuniform segregation of Impurities. 

These factors require that the furnace exhibit no radial temperature 
gradients and also that the G/R ratio be high. Finally, In the verti- 
cal Brldgman-Stockbarger technique the exact shape of the bottom of 
the ampoule Is important for the selection of a single seed from the 
many that may nucleate spontaneously. 

2. X-Ray Analysis 

Another important experimental technique employed in this investi- 
gation, has to do with the quantitative x-ray analysis of the Pb^ Sn Te 
rod grown by the Bridgman technique. Even before 1920, it was well known 
chat each chemical element emitted characteristic x-ray wavelengths when 
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( 24 ) 

excited by electrons or primary x-rays of sufficient energy. It was 

also observed that the intensity of the characteristic lines of an 
element depended on the amount of the element present in the specimen. 

So quantitative analysis by electron probe requires the measurement of 
wavelength and intensity of characteristic x-rays generated in the 
specimen. There are two systems used to separate the characteristic 
wavelengths for measurement. These are the wavelength dispersive 
spectrometer (WDS) and the energy dispersive spectrometer (EDS). Most 
of the analysis in this work was done with the energy dispersive 
spectrometer using a solid-state detector. The x-ray signal from the 
sample passes through a beryllium window and hits a Si (LI) detector 
generating electron-hole pairs. These are converted to voltage pulses 
and a multichannel analyzer sorts the pulses by voltages in proportion 
to their incident energies. In the wavelength dispersive spectrometer, 
an analyzing crystal is used to diffract x-rays from the sample in 
accordance with the Bragg's law. So at any given angle an x-ray of a 
particular wavelength is diffracted by the crystal and detected by a 
proportional counter. The x-ray quantum causes a number of ionizations 
in the detector which in turn creates a single output pulse whose 
amplitude is proportional to the number of ionizations. The total 
number of x-ray quanta is a measure of the relative number of atoms of 

the element giving off the particular characteristic x-ray quantum. 

( 25 ) 

It wasCastaing who first stated that the ratio of the characteristic 
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x-rays generated from element A In a specimen to that from pure A 
was equivalent to the concentration of element A In the specl'^z^n. 
That Is 


- <=A 


where ■ the Intensity of pure A x-rays 

s 

" Intensity of x-rays from element A In the sample for analysis 
" the concentration of A In the sample. 

However, this value commonly referred to as K value Is adjusted 

for several effects. These Include (1) the differences In electron 
scattering and retardation In the specimen and the standard; (2) the 
absorption of x-rays within the specimen, and (3) the x-ray fluorescence 
due to the different elements. This gives the following 
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( 22 ) 


where Kj,, K^, Kp represent the correction factors for the above mentioned 
effects respectively. 

C. EXPERIMENTAL PROCEDURE 

1. Furnace Characterization 


The furnace constructed for this investigation was characterized 
experimentally before it was used for the solidification experiments. In 
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order to determine both the axial and radial temperature profiles of 
the furnace, the temperature profiles in the center of the furnace, the 
wall of the furnace quartz lining, and the wall of the isothermal liner 
were measured. This was done using three thermocouple junctions lined 
up with a light microscope so that In a vertical position the junctions 
maintain equal furnace axial position but different radial position 
according to the furnace positions mentioned above. These thermocouples 
are run down the furnace starting from a recorded position with respect 
to the furnace reference point (FRF) using a motor and gear system. 

As the thermocouples are run through the furnace, a programmable Fluke 
thermometer records the temper^'ture at set time intervals. This -empty 
furnace characterization Is done first with both top and bottom ends of 
the vertical furnace kept open. It is then repeated with both ends 
plugged up with glass wool and some flberfax. Both the unplugged and 
plugged modes of the above experiment are done first without a tickler 
heater and then with each of the three tickler heaters described in the 
last section. A tickler heater was physically in the furnace at all 
times. So an ecperlment without a tickler heater simply means that the 
tickler heater power supply is cut off. It needs to be pointed out that 
the furnace design did not allow a complete temperature profile to be 
obtained between the Isothermal liner and the quartz liner (see Fig. 3). 
All three thermocouples for empty furnace temperature measurement are 
run through the hot zone. The thermocouple between the isothermal liner 
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and the quartz tube (liner) Is removed before the other two are continued 
Into the cold zone. The hot and cold zones heater controllers are set 
at 11.32 millivolts and 3.13 millivolts, respectively. This gives a 
temperature of 1100° C and 350° C respectively. These settings remain 
unchanged for all the experiments done with the furnace of figure 3. 

2. Melt-Solid Temperature Measurement 

The furnace was used to grow rods of Pbj^_^Sn^Te, Ge and Ag. The 
silver that was used In this Investigation was 99.99 pure but the 
other two matei.lals were 99.9999 pure. These were obtained commercially. 

A carefully weighed out portion of any of the three materials mentioned 
above is loaded Into an ampoule cleaned as described above and then 
attached to one of the vacuum pots. This loaded ampoule is pumped on 
as shown in figure 6, for about 12 hours. The weighing of the charges 
Is especially important in the case of the alloy semiconductor, Pb. Sn Te. 
The reason for this will be discussed in the next section of this chapter. 
However In all three materials, enough charges are weighed out to give 
about 18 cm of melt column in the 11 millimeter (ID) ampoule. After the 
loaded ampoule has been vacuum pumped for about 12 hours, it is then 

sealed under the vacuum. The vacuum gauge would typically read about 

•8 “6 
10 torr before sealing and about 10 torr just before the completion 

of sealing. A platlnum-platinum/rhodlum thermocouple threaded through 

a two bore alumina rod is slid down the ampoule capillary tube which has 

a 2-mm ID. This is shown in figure 5. The loaded ampoule, with the TC 
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in Che capillary, is lowered into the nearly isothermal region of the 
hot zone c£ the furnace as determined by Che empty furnace character- 
ization. This region is between 1100^ C to 1050° C. With a scale 
mounted on the furnace supporting beam, the position of Che bottom of 
the ampoule and the capillary TC junction are recorded with respect to 
the FRF. The whole system is then allowed to sit for at least 12 hours 
Co allow for melting, complete mixing in case of the Pb^^^Sn^Te and 
thermal equilibration. The TC in the ampoule capillary tube (henceforth 
..-eferred Co as capillary TC) is connected to a programmable thermometer 
by Fluke. This Fluke unit Is then programmed to take temperature readings 
at some time interval— usually ranging from 30 seconds to 30 minutes. 

The ampoule and the capillary TC are lowered towards the cold zone by 
a drive motor with a predetermined velocity. The capillary TC measures 
Che temperature in the melt at the set time intervals as the system is 
lowered. Eventually, a solid phase forms and the advancing solid-melt 
interface in due time reaches the capillary TC junction. As Che ampoule 
is lowered further, the capillary TC starts to measure the temperature in 
the solid phase at the same set time intervals. The temperature measure- 
ment usually continues for about 1 centimeter into the solid phase as 
determined by the distance scale on the furnace. Then Che capillary TC 
junction is pulled up manually back into the melt as determined by the 
temperature reading. This is done while Che ampoule is still being lowered. 
So at no time is the anpoule movement stopped. On completion 'f the 
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capillary TC pull-up the thermocouple Is clamped down and allowed to be 
lowered again together with the ampoule. This makes It possible to 
obtain another set of temperature measurements. Each time this experi- 
ment Is repeated, the thermocouple reads the temperature profile from 
the melt phase through the melt-soHd Interface Into, the solid phase. 
Typically this p rocess is repeated at least five times for any run. The 
experiments described above in which the ampoule Is In continuous motion 
are referred to as dynamic runs. 

On completion of the above described dynamic runs, the drive motor 
Is turned off. The system Is allowed to equilibrate for at least 4 hours. 
This allows for the dissipation of heat generated due to the latent heat 
of fusion and/or heat carried due to ampoule movement (Peclet effect) . 

A constant temperature reading over a period of about 30 minutes was 
taken as Indication of complete thermal equilibration. With the ampoule 
held stationary, the temperature measurement experiment is repeated by 
running the capillary TC down the capillary tube. The tickler heater 
setting is lowered and the system is again left to equilibrate before 
another temperature measurement. All the temperature measurement experi- 
ments done with the ampoule at one position are referred to as stationary 
runs. The important variable in the stationary runs Is the tickler heater 
millivolt power supply setting. A example of Pb^^ ^Sn^Te and Ag direction- 
ally solidified as described above is shown in figure 7. 
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Figure 7: Directionally Solidified Rods of Pb^ 
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3. Quanfltatlve X-Ray Analysis 

The alloy semiconductor, Pb, Sn Te, grown by the above procedure, 
was analyzed for compositional variations along the grown crystal rod. 

This was done using quantitative energy dispersive x-ray analysis. The 
energy dispersive analysis was done using a JOEL c'*a'nlng electron 
microscope (SEM) in conjunction with a TRACOR TN 2000. Binary alloy 
standards FbTe and SnTe were used to collect reference spectra at 
recorded beam current, count rate and total count time. The electron 
beam current used in all the analysis is set using a Faraday cup for 
beam collection and an ammeter for signal readout. The x-ray spectra 
from the samples taken from the grown rod of Pb._ Sn Te are collected 
under the same conditions mentioned above. Other conditions that 
remained the same in all the analysis are the take-off angle and detector 
to sample distance. These x-ray spectra are then analyzed for compositions 
using the "super ML" program of the TRACOR unit. The unit alro has a 
"ZAP" program for the corrections for atomic absorption or enhancement, 
matrix and fluorescence effects. 

D. SAMPLE PREPARATION 

1. Compounding Pb.,_ Sn Te 

Lead and tin combine with tellerium to form binary compounds PbTe 

( 26) 

and SnTe. These possess the cubic Bl (rocksalt) crystal structures. 

The composition range within which PbTe and SnTf* and exist, as showi, in 
figure 8a, i? so narrow that they can be represented as compounds with a 
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single line at 50 atomic percent. Also shown in figure 8a is the 
fact that there is some range to the composition within which the 
compound can exist. In this investigation the ternary systems Pbj^_^Sn^Te 
is treated as a pseudobinary system. as shown in figure 8b. This is simply 
the shaded isopleth of figure 8c and is based on the assumption of a 
50 atomic precent pure compound PbTe and SnTe. It was therefore very 
important to weigh out the elements carefully enough to have a one to 
one ratio between the metals (Pb,Sn) and nonmetal (Te) in order to maintain 
stoichiometry. Once stoichiometric proportions are maintained it becomes 
a good approximation to treat the system as an isomorphous pseudobinary. 

The appropriate stoichiometric weights of the elements are loaded in a 
quartz ampoule, melted in a vertical furnace and left to react for about 
12 hours before the solidification experiment. Sometimes the charge is 
melted in a rocking-furnace which may enhance proper mixing and re<i..tion. 

2. Cleaning and Etching of Samples 

Both the germanium and the silver charges used in this work were 
cleaned by a chemical etc. The cleaning solution for the germanium was 
as follows: 

3:5:3 by volume of HF: HNO^t Acetic acid plus a drop of Bromine. 

The germanium charge was left in this solution for 2 minutes and 

then rinsed copiously with deionized water. The silver charge was 

cleaned in a solution of ammonium hydroxide (NH, OH) and hydrogen 

4 

peroxide (H 2 O 2 ) in a 1:1 ratio by volume, using the standard laboratory 
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stocks. The cleaning takes about 1 minute after which the silver is 
rinsed in water and then in acetone. Typically, the quartz ampoules 
in which these materials are grown, are cleaned as follows: 

a) Clean with trichloroethlene 

b) Rinse with methanol 

c) Rinse with deionized water copiously 

d) Etch for 10 minutes in a mixture of 60S deionized water, 

35Z nitric acid and SS hydrofluoric acid all by volume 

e) Rinse with methanol copiously 

f) Torch-dry with oxy-acetylene flame 

g) 1000° C bake-out undei a vacuum system 

A grain boundary etch was done on the grown Pb^ ^Sn^Te crystal using 
the Norr solution.'*' The composition ot the etchant is 20 gm of KOH 
in 45 ml of H 2 O; 35 ml of and 20 ml of C^H^OH. This electrolytic 

etch is done with 1 volt potential for 15 secs and then rinsed 3 to 4 
times in ethanol. 

3 . Sa mples for X-Ray Analysis 

In Che preparation of samples for the quantitative x-ray analysis, 

both mechanical and chemical polishing was used. The Pb, Sn Te rod is 

1-x X 

carefully sectioned into samples for analysis using a diamond wheel cutter. 
These samples are polished mechanically starting with a dry 400 grid 
emery paper and then taken progressively through the finer papers. The 
final mechanical polish is done on clothe mounted on steel wheels using a 
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1 micron alumina abrasive first and then a 0.05 micron alumina abrasive. 
It is possible to get smearing of the material on the polished surface. 
Tj eliminate this possible smearing problem, the mechanically polished 
samples are given a 1 to 2 minute chemical polish. The chemical polish 
is 95Z HBr and 5t Br by volume. 
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CHAPTER IV 

EXPERIMENTAL DATA AND RESULTS 

A. Introduction of System Geometry 

In order to follow the data and results to be presented In this 
chapter, it is important to introduce the system geometry and defini- 
tion of coordinates adopted in this investigation. The system 
geometry and coordinates are shown in figure 9. A line half way along 
the insulation zone is taken as the furnace reference point (FRP) , 
from which all furnace axial positions are measured. The symbol Z 
represents any position along the furnace axis. The symbol i denotes 
the position of the melt-solid interface with respect to the furnace 
reference point. Both Z and i have a zero value at the FRP, are 
positive toward the cold zone and negative toward the hot zone. If 
the bottom of the ampoule is located at a particular value of Z, 
say Za, then the length, L, of solidified material is given by 

L - Za - i (23) 

and the fraction of the melt column solidified at any given instant, 
gg, is to first order (Za - i)/L^ where L^ is the total length of 
the boule when it is completely liquid. 

B. Furnace Characterization 

The empty furnace characterization vas accomplished by pulling a 
set of thermocouple junctions through the furnace of a known pull rate 
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while a digital thermometer prints the temperature of the thermocouples 
at a preprogramned time interval. The raw data obtained from the 
digital thermometer is in the form of temperature versus time. These 
data are fed into a computer which converts time to distance, z, and then 
generates temperature and/or temperature gradient versus distance plots 
by use of the programs shown in Appendix I. 

The empty furnace temperature profiles are shown in figures 10 and 
11 for the case of the tickler heater's power supply cutoff. That means 
that the tickler heater was in position in the furnace but was not 
generating any heat. The difference between figures 10 and 11 is that 
for figure 10, both the hot and cold zone ends were open to air flow, 
whereas both ends were plugged with quartz wool in the case of figure 11. 
Notice in figure 10 the large difference between the near-wall and the 
center-line temperatures, up to 300° C. Figure 11 shows that plugging 
the furnace ends results in a remarkable reduction in the temperature 
dilferences between the near-wall and the center line, less than 100° C. 

Nc tice also that the maximum differences in near-wall and center-line 
temperature occurs at the boundary of the hot zone and the insulation 
zone. The above mentioned figures also show that the center- line tempera- 
ture gradient increased from 22° C/cm at 800° C to 75° C/cm at 800° C 
in the plugged furnace. The furnace design did not permit the temperature 
measurement near the wall of the isothermal liner to continue through to 
the cold zone. So only the hot zone profiles are shown in the figures 
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EMPTY FURNACE TEMP. PROFILE 



Figure 10: Unplugged temperature profile with a 2.3 cm tickler 

heater off. A few minutes after plugging the curve 
changed to that in Fig. 11. 
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mentioned above. If a possible furnace misalignment Is neglected and 
a perfectly synnetrlc furnace geometry Is assumed, figure 11 can be 
used for the construction of Isotherms along the furnace as Is shown 
In figure 12. This figure Is typical of most of the empty furnace 
temperature profiles In the sense that the Isotherms remained convex 
(looking from the hot zone) all through the furnace but with a decrease 
In convexity as It approaches the cold zone. As will be noted later 
this Is not the case for a loaded furnace. 

Measurements of the plugged^fumace temperature profile was 
repeated but this time with a 1.43'-cm long tl.ikier heater activated to 
hold a reference thermocouple at about 1045° C. The resulting tempera- 
ture profile is shown in figure 13. This figure shows a center-line 
temperature gradient of 87° C/cm at 800° C which means an Increase in 
temperature gradient over that of figure 11. Figure 14 Is a plot of 
the temperature gradient exhibited by the quartz liner's near- wall 
profile shown in figure 13. It shows a near -wall temperature gradient 
of 120° C/cm at 800° C as opposed to the 87° C/cm mentioned above for 
the center-line. 

A much better plugging of the furnace ends was later achieved by 
using quartz wool first and then covering it with fiberfax (alumina) 
specially grooved to fit the furnace ends better. This relatively 
better plugging resulted in the temperature profile of figure 15 which 
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' IV- 4 

shows an almost complete elimination of temperature differences between 
the near-wall and the center- line profiles. This means that any radial 
temperature gradient is almost removed by effective plugging. Figure 15 
also shows that the axial temperature gradient was increased to 130° C/cm 
due to a combined effect of a 'longer tickler heater (2.3 a' and a 
relatively good plugging. 

It is known that the presence of a load in a furnace shifts the 

2 

isotherms which means a shift also in the temperature profile. The 
conduction of heat by the load tends to flatten the furnace profile. 

The motion of the load pushes heat toward the cold zone and therefore 
shifts the isotherms. The temperature difference between the near- 
wall and the center-line of a loaded furnace as shown in figures 16 
and 17a Is very small when compared with that of the \mplugged furnace 
profiles of figure 10. It is noticed that in figxires 16 and 17a, for 
a given position the center-line temperature is lower than the near-wall 
temperature up to a point. After this point the center-line temperature 
becomes higher and remains higher than the near-wail temperature for the 
rest of the furnace. Of course, this is in contrast with the empty 
futnace profiles in which the center-line temperature remained lower than 
the near-wall temperature all through the furnace. Employing the same 
assumptions made for figure 12, the loaded furnace profile can be con- 
verted into isotherms as shown in figure 17b. However, in contrast with 
figure 12, the isotherms change from convex (looking from the hot zone) 
to concave, going toward the cold zone. 
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♦Furnace reference point. 



63 


C. Melt-Solid Temperature Measurements 

1. Interface position determination 

The raw data from the meit-aolid experiments described in Chapter III 
is in the form of temperature versus time. These data are reduced to 
temperature versus distance and temperature gradient versus distance 
along the furnace. Curves based on these reduced data show a discon- 

4 

tinuity in slope at the melt-solid interface. For each of the dynamic 
or stationary experiments described in Chapter III, a plot of temperature 
and/or temperature gradient against position, z, in the furnace is made 
to determine the interface position. It is obvious then that many curves 
of this type were in this investigation, however, in order to avoid clutter, 
only a few examples of each type will be given here and in Appendix II in 
order to represent the three materials used in this investigation. 

The complete set of data points obtained in one of the dynamic 
experiments is shown. in figures 18 and 19. Both figures portray a discon- 
tinuity in the regions marked AB. However, the derivative of Che tempera- 
ture with respect to distance is a more sensitive condition as indicated 
bv the more pronounced discontinuity in figure 19. On expanding figure 18 
within the region marked /\B, the interface position with respect to the 
furnace, can be determined more clearly as shown in figure 20. This gives 
Che first interface position monitored for a boule of Pbj^_^Sn^Te and for 
the indicated conditions. The freezing point of the material is deter- 
mined bv noting the temper.sture at the point of discontinuity. The curves 
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Figure 18; Dynamic temperature measurement from the melt into 
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for the second through the fifth interface positions monitored in this 
particular run are in Appendix II. Also shown in this Appendix are 
examples of the same type of curves generated for the interface position 
and freezing point determination for silver and germanium. 

2. Materials constants 

From the experimental technique employed in this investigation, two 
material constants can be determined. These are the freezing point anJ 
the ratio of the melt to solid thermal conductivities at the freezing 
point. It has been mentioned that a plot of temperature versus distance 
from the melt into the solid, exhibits a discontinuity in slope at the 
melt-solid interface and that the temperature corresponding to this 
point of discontinuity is the freezing temperature of the material. 

From figures like figure 20, the freezing poii ts of silver and germanium 
were determined from an average of about 15 runs each to be 960.1 +0.9° C 
and 937.3 + 1.5° C respectively. The uncertainty values are one standard 
deviation from the mean. These values are in good agreement with the 
literature values of 961.9 C and 937.4° C for the respective freezing 
points. As would be expected for ^n alloy material which changes com- 
position during freezing, the freezing point of Pb^ ^Sn^Te varied as 
freezing progressed. The measured freezing points and the corresponding 

compositions will be indicated later in this section. 

2 3 

Theoretically, ’ it is known that the thermal conductivities in 
the solid and melt phases influence the temperature profile during a 
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directional solidification (DS) experiment. Table 1 shows the measured 
temperature gradients at 900*C for the three materials, lead-tin- 
tellurlde, germanium and silver. Note that for the. same furnace tem- 
perature settings, the temperature gradient is reduced as the thermal 
conductivity of the load material Increases. The thermal conductivities 
of the solid phase, Kg, Of many. materials are easy to find in the litera- 
ture whereas those of the melt phases, especially of high melt-point 
materials, are very dlffucult to find. However, proper heat balance ana- 
lysis across the Interface in a DS experiment can in theory yield the 
thermal conductivity ratio of the melt and solid phases as the freezing 
point. Assuming a planar Interface and ignoring both radiative and con- 
vective heat transfer, the equation for the heat balance Is 


+ pRH - KgGg (24a) 

where p ■ density 

R “ melt-solid interface velocity 

H » latent heat of freezing 

S " temperature gradient in the liquid and solid, respectively 

Furthermore, if a stationary experiment is considered, that is R “ 0, 
then the equation reduces to 








(24b) 


Since Gg/G^^ 

determined. 

32 

of silver 


can be determined from figure 20, can also be 

Thermal conductivities of both the solid and liquid phases 
33 

and germanium are recorded in the literature. So part of 



69 


Furnace Condition 

Temperature 
Gradient at 
goo^’C 

Tickler 
Heater 
Setting 

Empty 


127°C 

1045“C 

Loaded 
Wi th 

PbSnTe 

118“C/cm 

1045°C 

Ge 

gg^c/cm 

1045°C 


Ag 

3G®C/cm 

104b'’C 


Table 1: 


Effect of a load material on the furnace temperatu'^e 
gradient at 900°C. A 2.3 cm tickler heater was used 
in all the cases. 
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the reason tor performing the DS experiments with silver and germanium 
was to check the accuracy of the experimental technique in determining 
the ratios. Table 2 shows the measured 

K^/K^ ratios as compared with selected literature values. The com- 
parison indicates that the measured value is too low for both the 
germanium and the silver. The error in the thermal conductivity ratios 
will be discussed semiquantitatively in the next chapter. In any case, 
the ratio measured for Pbj^ ^Sn^Te is 1.75 + .04. 

3. Tickler heater effects 

A tickler heater (booster heater) wao introduced into the furnace 

in order to increase the melt temperature gradient which is known to 

4 11 

affect interface stability. ’ The extra heat generated near the melt- 
solid interface by the tickler heat can also influence the position of 
the interface with respect to the furnace. These possible effects were 
studied by currying out the melt-solid temperature measurement experi- 
ments with each of three tickler heaters of different lengths. 

The interface positions determined from any of these experiments 
were used to plot a graph of interface positions versus the translated 
length as shown in figure 21- 

The figure shows that for the two materials Ge and Pbj^_^Sn^Te, 
the melt-solid interface position did not significantly move with 
respect to the fixed furnace reference point within the solidification 
period that was monitored. For the same material and the same tickler 
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Figure 21 : 


Crystals grown ith the same ampoule translation 
rate but with different tickler heater lengths 
and settings. 
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heater power supply setting the interface position was depressed toward 
the cold zone as the tickler heater length increased from 1,43 cm to 
2.54 cm. The interface position is least depressed when the tickler 
heater is cut off. The germanium melt-solid interface is more depressed 
than that of Pb, Sn Te. Both were solidified using a 2.54 cm tickler 
heater but different power supply settings. It was set at 1095*^ C for 
germanium and 1045° C for Pb^^ ^Sn^Te. All the interface positions 
shown in the figure were determined from dynamic runs, and the tickler 
heater setting was constant all through each run. In any case, it can 
be seen from the above results that the shift in interface position can 
be due to the tickler neater length. The tickler heater setting can 
alsc shift the interface position. A direct demonstration of interface 
position shifts due to tickler heater temperature setting is shown in 
figure 22. These interface positions are determined from stationary 
runs. The only cause of solidification in this case is the change in 
the tickler beater power supply setting. This shows that for a given 
material tne interface position shifts towards the hot zone as the 
tickler heater setting temperature is lowered. Within the experimental 
error, the interface position varied almost linearl” with the tickler 
heater temperature setting. 

Temperature gradient in the melt Immediately ahead of the interface 
is an important factor in crystal growth from the melt in general and 
especially in the Bridgman-Stockbarger technique. This is particularly 
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EFFECT OF TICKLER-HERTER SETTING 



TICKLER-HEHTER SETTING (“O 


Figure 22: Effect of tickler setting on interface position. This is 

a stationary run so all interface movement was due to a 
reduced tickler heater setting. 
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true in the growth of alloys which require the avoidance of both thermal 
and constitutional supercooling. . Table 3 shows an increase in melt 
temperature gradient with an increase in tickler heater length for the 
furnace configuration employed in this investigation. For the seme 
tickler heater length and temperature setting the temperature gradient 
in the melt is lowest for silver which has a relatively high value. 

In conclusion, for the furnace constructed and used in this investi- 
gation, the temperature gradient in the melt immediately ahead of the 
interface can be adjusted by a change in the length and/or temperature 
setting of the tickler heater. These changes also result in shifts in 
the interface position. 

A. Ampoula translation rate effect 

Two translation rates, .0A6 cm/min and .178 cm/min were used in the 
DS experiments. The effect of the ampoule translation rate V on the 
melt-solid interface position of a material, viz Pb^ ^Sn^Te, is shown in 
figure 23. The interface shifted toward the cold zone with an increase 
in translation rate. This effect is expected to be true also for both 
silver and germanium. However, the degree of the shift may be different 
for different materials. This will be explained in the next chapter. 

A very important result shown alsf) by figure 23 is the fact that as 
more material is solidified (higher L value), the melt-solid interface 
position, i, is constant for tb«* sem:. conductor , Pb._ Sn Te. This is also 
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Material 

Tickler Heater Length* 
and 

Power Supply Settings 

Temperature 
Gradient in 
the Melt Ahead 
of the Inter- 
face 

PbSnTe 

2.3 on- 
cut off 

72* C/cm 

1 .43 on 
1045“C 

80* C/on 

2.54 cm 
1045°C 

37° C/cm 

Ge 

2.3 cm 
1045°C 

50* C/cm 

2.54 cm 

1095°C 

63° C/cm 

Ag 

2.3 cm 
1045'’C 

35* C/cm 

2.54 cm 
1045®C 

40° C/cm 


*6 heater wire windings per cm. 


Table 2: Temperature gradient in the melt immediately ahead of the 

interface showing an increase with tickler heater length. 
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Figure 23: Rods solidified at tw^ Afferent ampoule translation 

rates but vAth a 1.43 l..< tickler heater set at 1045°C. 
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true for Ge as was already mentioned in connection with figure 21. 
However, for a metal, Ag, there is axi appreciable variation in the 
interface position as freezing progresses. In fact, it moves toward 
the hot zone (higher negative i values) and apparently in a nonlinear 
manner. This is interpreted to mean a shift in the isotherm that 
corresponds to the freezing temperature which in turn means a continu- 
ous change in the temperature profile that the material sees. The 
different temperature profiles seen by four interfaces monitored in 
one DS experiment with Ag are shown in figure 24. The profile shifts 
toward the hot zone as freezing progresses and hence a decrease in the 
temperature gradient in the melt at the freezing temperature also occurs. 

5. Growth rate variations 

The metal, silver, showed a continuous variation in tb® instantan- 
eous growth rate as freezing progressed. The semiconductors, Pbj^_^Sn^Te 
and Ge shc:,ed a constant growth rate within the growth period monitored. 
These results are deduced from figure 25, since the frozen length, L, 
is directly dependent on the translated length Z. The crystal growth 
rate R ■ dL/dt where dL/dt is given by differentiating equation (23) 
with respect to time, t 


R 


dZ 
a 

dt 


di 

dZ 

a 


dZ 
a 

dt 


V(1 - f) 


( 25 ) 


where the ampoule translation rate V ■ dZ /dt and dZ - dZ. 

3 3 
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This equation simply means that the growth rate Is the sum of the ampoule 
translation rate, V, and t’ amount, dl/dt by which the growth rate 
differs form the translation rate. Note that the difference between the 
Instantaneous growth rate R and the ampoule translation rate, V which 
is proportionate to the derivative of the interface position with respect 

to the translated length, was close to zero for both Ge and Pb. Sn Te. 

X*X X 

However, for silver this difference was initially high; then it became 
progressively smaller and at a critical point, it turned up and became 
progressively higher. Using equation (25) it can be seen from figure 25 
that the growth rate of silver remained higher than the ampoule trans- 
lation rate all through the growth period shown. At one point the silver 
growth rate was higher than the translation rate by as much as 110%. The 
growth rate of Ge and Pb^ ^Sn^Te was equal to the translation rate all 
through the growth period monitored. 

It is expected from equation (25) that if the crystal growth rate, 

R, and ampoule velocity are constant, the fraction of material solidified, 
L/L^ will have a linear dependence on the translated length. The fraction 
solidified is the length solidified normalized with the total boule length. 
The translated length, Z, Is made dimensionless by dividing with the 
diameter, d, of the boule. Figure 26 shows a plot of L/L^ against 
Z/d or each material tested. For the semiconductors, Ge and Pb, Sn Te , 
the relationship is essentially linear within the growth period shown. 

For the metal, silver, the relationship is nonlinear as could be expected 
from the curve fcr silver in figure 25. 


di/dz 
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Figure 25: Growth rate variations due to the length of material 

sol idified. 
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D. Composition Analysis 

As already mentioned in this chapter, the freezing temperature for 
the alloy semiconductor, Pbj^_^Sn^Te, varied as more of the material was 
frozen. This is in agreement with what would be expected for a pseudo- 
binary Isomorphous system. A qualitative x-ray analysis was performed 
on the grown rod crystal to determine the solute profile and consequently 
to determine the compositions that yielded the various freezing temperatures 
measured in this investigation. 

Figure 27 shows the concentration profile of the crystal rods grown 
at the indicated translation rates. The curve is fitted to nquation (12) 
in Chapter II, but with the effective segregation coefficient, sub- 

stituted for the equilibrium segregation coefficient, k, which gives 


C 


-1 


k C 
eff o 


(1 - 83) 


eff 


( 26 ) 


This is done because equation (12) only applies to the particular case 
where there is complete mixing in the liquid and no diffusion boundary 
layer. However, if the k, in equation (12) is changed to then 

the equation becomes a more general one for a system with mixing in the 
melt. Note that from the definition of equation (10) the general 

equation reduces to the particular equation (12) when 6 is zero where 
6 is the thickness of the diffusion boundary layer. From the general 
equation (26), the ^gff calculated for both composition-distance 
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Figure 27; Concentration profiles of Pb, ^Sn Te crystal rod grown 
at different translation rates, it also shows the mea- 
sured solidus temperatures. -1 

■^Calculated from the equation C = *^eff^o^^ " 




85 


curves, giving a value of .790 for a translation rate of .046 ca/mln and 
.810 for a translation rate of .173 cui/mln. These values are obtained 
by using the above equation (26) In the logarithmic form namely, 

InC • Ink^j^Co + - 1) In (1 - g^) (27) 

The least squares fit to the above equation (27) for the entire data gives 
a straight line whose intercept is Ink^^^Co and slope is (k^££ •" D* From 
this, hot!) Co and ^^££ can be calculated. The Co calculated by this 
method is .155 mole fraction SnTe for the material grown at .046 cm/min 
and .215 mole fraction of SnTe for that at .178 cm/mln. These values 
compare favorably with the value 0.2 mole fraction determined by weighing 
during sample preparation. 

Also shown in figure 27 are some of the interfacial temperatures 
measured at the same Z values for which >"he composition measurements 
were made. This temperature composition data in effect represents an 
effective nonequilibrium phase diagram which will be discussed in 
Chapter VI. Notice here that the temperature composition data closely 
correspond to those numbers obtained from the liquids line of the pseudo- 
binary equilibrium diagram of figure 7. 

A typical x-ray spectrum collected and used for the quantitative 
analysis is shown in figure 28. It gives the M and L x-ray lines 
of lead and L lints of both tin and tellerium. On carrying out a grain 
boundary ^'■('h on this sample using the method described in Chapter III, 
a jiicrography of the surface reveals some grain boundaries as shown in 
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Figure ?8; tnergy dispersive x-ray spectrum cf a cross-section of a 

Pb, Sn Tc rod. This is an example of the spectra used in 
th^’^uat^ti tati ve x-ray analyses. 
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figure 29. Also shown in figure 29 is an x-ray spectrum typical of the 
grain boundaries. The Fb and Sn lines are almost completely absent from 
the grain boundary x-ray spectrum. This may mean that the grain boundary 
etch preferentially etches out the lead and tin leaving tellerium etch 
botmdarles. No similar spectrum of mostly tellerium was observed in 
the unetched specimens. 

A sample of the EDS results plotted in figure 27 for V » .046 cm/min 
is shown in Table 4 in terms of the atomic percentages of the elements. 

It is noted that at every point analysed, the sum of the atomic percentages 
of the metals, Pb and Sn, is always approximately equal to the atomic per- 
centage of the nonmetal, tellerium. This may Indicate that the crystal 
grows by the incorporation of stoichiometric molecules of PbTe and SnTe 
as opposed to the atcxos of the elements. 

In line with the above resiilt, it was important to ensure that grains 
of pure lead, tin or tellerium were not included in the grown sample due 
to poor mixing and/or non-reacLion of the elements. As already described 
in Chapter III, the elements are weighed out in a proportion that should 
give stoichiometric Pb^ ^Sn^Te when properly mixed and reacted for about 
twelve hours in the furnace. Several samples were taken from the grown - 
material and powdered for x-ray diffraction experiments. Table 5 shows 
one result of such an experiment. The data shows no d spacing tiat matches 
any of the three elements, as seen by comparison with the data fcr the 
pure elements taken from the standard ASTM x-ray index cards shown in 



Table 6. 


Position Along the 
Crystal Rod {cm) 


5.3 


6 . 


19.13 


10.52 


12.94 


14.33 


Atomic Fraction of the 

Elements 

Pb 

Sn 

Te 

.398 

.093 

.508 

.409 

.084 

.507 

.405 

.088 

.507 

.397 

.097 

.506 

.385 

.106 

.508 

.394 

.101 

.505 


Table 4: Result of quantitative x-ray analysis of the interface 

positions monitored during the Pb, Sn Te growth at 
V = 0.46 oi/min. 1-x x ^ 

Note that the composition is always approximately .5 
metal (Pb and An) a'.d .5 nonmetal (Te). 
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OF POOR QUALITY 

From Table 5» the lattice constant of the sample can be calcualted 
using the formula for a cubic system and was found to be .64 nm. 

It Is assumed, at least for the first to freeze pdrtlon of the 
crystal rod (low sclute concentration), that Vegard's law Is obeyed. 
Based on the assumption, the composition of the sample can be calculated 
as shown In Appjnd^ I.* . results oi this calculation Is .143 mole 
fraction of Bnte whlcb-compares. favorably to the composition measured by 
the x-ray analysis result shown In figure 27 for the first to freeze 


region. 
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Figure 29: The x-ray spectrum shown is typical of as etched arain 

b&undar i -■'S of Pb, Sn Te shown in the upper part of the 
figure, ^he grain boundaries are depleted of Pb and Sn. 
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CHAPTER V 
ERROR ANALYSIS 

A. Melt-Solid Thermal Conductivity Ratio 

The thermal conductivity ratios are determined from a measurement 
of the temperature gradients present In the material In the growth fur- 
nace. The errors In conductivity, therefore depend on the errors In 
the thermal gradients which can be assessed from the scatter In the data. 
The temperature gradients of the melt and solid respectively, and Gg 
used In the calculation of the thermal conductivity do show a scatter In 
the data. This scatter reflects the time dependent error In the position 
reading and the error In the temperature reading. There Is no easy way 
to assign an error to the temperature reading. However, given the fact 
that the freezing point measured for the pure samples, Ge and Ag, were In 
good agreement with the literature values. It can be assun^d that the 
error In the temperature reading Is small. The scatter In G^ and Gg 
measured at points near the Interface Is taken as the error In the tem- 
perature gradient measurement. This ranged from 2 to 5Z and therefore 
cannot explain the 33Z average error reported In Table 2 for the deter- 
mined thermal conductivity ratios. This discrepancy Is Interpreted to 
mean that systematic error Is present and several possible sources of 
systematic error will be discussed below. 

It was mentioned In the chapter on experimental details that the 
thermocouples had to be protected from chemical attack during the 
melt-solid temperature measurement. This was achieved by running the 
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thermocouples through a qxuirtz capillary tube. Tlie fact that this capil- 
lary tube runs through both the melt and solir’ phases will Introduce 
a systematic error in the thermal conductivity ratio measured for the 
materials. The capillary tube a;ts as a conductor in parallel with the 
material to be solidified. Therefore, the calculations of thermal con- 
ductivity from equation (24b) which utilizes temperature gradient mea- 
surements in actuality yields the thermal conductivity ratios of the 
glass and sample rather chan that for the sample alone. The magnitude 
of error in the thermal conductivity determined in this manner due to the 
presence of the capillary tube is considered using figure 30. 

For the tube plnp molten material, one has a net conductivity K 

tni. 


where iL * thermal conductivity of melt, glass 

F - volume fraction of melt, glass, 

with F, + F ■ 1. 

L g 

Similarly for Che solid side of Che interface, the net conductivity 


+ Vs 


(29) 


where K- and F are Che thermal conductivity and volume fraction respec- 
d S 

tively of the solid. From the dimensions of the ampoule (Figure 5) it 


can be seen that 
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OF POOR QUAL: i t 


- Fg - .6214 (30a) 

F - .3786 (30b) 

S 

where the qoartz capillary and the thermocouple assembly was assigned the 
thermal conductivity of glass. The calculated thermal conductivity 
ratios f”om the measured temperature gradients are from equation (24b) 



K 


'mS 


G 


8 



(31) 


Substitution nf equations (28) (29) and (JOa) Into equation (31) 
yields 


Vs Yg =L 

which can be rearranged to yield 


■S. 


K 


b 


G_ K F 



(32a) 


(32b) 


The second term the right ..and side of equation (32b) Is the 
error Introduced by the presence of the capillary tube and can be ex- 
pressed as a fraction error by 


K T 
•nL 



K F G - 
_ 8 g L S 

K , 

riiO 


^mS 

‘=s’'s =s 


( 33 ) 
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For a given measurement of temperature gradient one can see from equa- 
tion (30) that the larger the product larger the fractional 

r. (32) , (33) . , (34) ^ ^ 

error. For silver , germanium , and gliss one has K_ ■ 

Ag 

3.55 watt/cm-°K, K- ■ 0.465 watt/cm-°K, K ■ 0.86 watt/cm-°K, 

Sge 8 

/Go - 1) - .355, and F„/F-(G. /G. - 1) - -.274 to yield 

* Ag Ag ® Ge Ge 

from equation (33) respective errors of about IX and -5%. 

The above calculated errors are very small compared to the 33% 
average error In the actual data shown In Table 2. It was, therefore, 
attempted to explain the error by considering possible geometrical ef- 
fects du£ to the capillary tube and furnace geometry. 

The effect of the possible curvature of Isotherms or the presence of 
a radial gradient was considered; figure 31 shows a possible arrangement 
of curved Isotherms that may exist In the loaded furnace. Note that the 
temperature measurements are made along the longitudinal axis which Is 
likely not normal to the Isotherms. However, equations (24a) and (24b) 
hold when the temperature gradient Is measured normal to the Isotherms. 
The normal temperature gradient G^^ is obtained by the vectorial addition 

of the measured axial gradient G„ and the measured radial gradient G . 

n r 

That Is 

G„ - /G^ + G^ (34) 

N m r ' ' 

From the temperature measurements along the ampoule wall and along the 
capillary average axial as well as average radial gradients can be 
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determined. The normal gradients for melt and solid and were 
calculated from equation (34) and substituted into equation (24b) to 
calculate ratios. This approach yielded thermal conductivity ratios 

which differed from the literature values by a larger amount than that 
for the uncorrected data in Table 2. 

All the above analyses failed to provide a suitable correction fac- 
tor that could be applied to the data. However, it is noted that the per- 
cent error between the values obtained here and those available in the 
literature was approximately the same for both germanium and silver as 
shown by Table 2. In view of the lack of a physically based correction 
procedure, this percent error (33%) was used as an empiiical correction 
factor and applied to the PbSnTe results to obtain what is likely a better 
assessment of the actual thermal conductivity ratio. This empirical ap- 
proach gave a Kj^/Kg value of 2.35 ± .08 for the Pbj^_^Sn^Te, boule where x 
is 0.2 for the case studied here. 

There is one additional comnent to be made about the use of equation 
(24b) to calculate the thermal conductivity ratio. When the boule and 
interface position are stationary, this equation is applicable; how- 
ever, when the boule is being translated at a velocity V and the solid is 
growing at a linear rate R, then the latent heat H may make a significant 
contribution to the heat balance equation. 

The heat energy balance equation at the interface is given by equa- 


tion (24a) as follows 
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+ pRH - KgGg <24a) 

whern p Is the density of the solid. Since the Gj^ and Gg are determined 
experimentax and the material constant are known, it la possible to sub- 
stitute these values Into the above equation in order to quantitatively 
determine the relative magnitude of the latent heat term In comparison 
with the heat conduction terms. According to the results, the growth 
rate of silver, R, is twice the ampoule translation rate. For a sample 

calculation, a translation rate of .178 cm nln"^ is used. This gives a 

growth late of 5.9 x 1C cm sec . 

*5 O 

For silver, H - 26.5 cal/gm, p - 10 gm/cm"’, • .4183 cal/cm-sec- C, 
K- - .8485 cal/cm-sec-°C. For one of the above dynamic experiments run 
at the above mentioned translation rate, the measured G^^ and Gg values 
are 40.44°C/cm and 25.56°C/cm respectively. Therefore, 

■ 40.44*^C/cm x .4183 cal/cm-sec-”c ■ 16.92 cal/cm^-sec 

K-G - 25.56°C/cm x .8485 cal/cm-sec-°C •• 21.86 cal/cm^-sec 

pRH » 10 cm/cv.^ X 10“^ cm/sec x 26.5 cal/gm " 1.57 cal/cm^-rec 

The total rate at which thermal energy flows into the Interface is 

iLG- + pRH ■ 18.49 cal/cm^-sec 


The latent heat of fusion term comprises only about 8.5Z of this heat. 
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For germanium, a similar calculation can be made using the measured 
and G. values of 58.74°C/cm and 107,38°C/cm respectively. The other 

Ij w 

relevant parameters for the germanium are p ■ 5.323 gm/cm , R ■ .046 cm/min, 

• .311 cal/cm-sec-°C, Kg “ .111 cal/cm-sec-°C. 

TMs gives 

KG - 18.25 cal/cm^-sec. 

L J ) 

C_K “ 11.93 cal/cm^-sec 
pPH ■ .466 cal/cm^-sec 

Therefore, the latent heat term in this case is about 2.5% of the total 
rate of heat flow into the n«lt-solid interface. Here also the latent 
heat term is considered negligible for the translation rates used in this 
work. 

However, notice from the above calculations that for silver 

K^C^ + PRH < K^Gg (35a) 

and for germanium 

Kj^Gj^ + pRH > KgGg (35b) 

These inequalities suggest that there is another heat transfer 
mechanism operating at the interface. Convection in the melt is probably 
the other heat transfer mechanism. In order to account for this effect, 


H - 114.3 cal/gm, K^^ 


one can consider an effective heat t ansfer coefficient for the melt, h , 



100 


which would include all the heat transfer coefficients regardless of Ae 

transfer mechanism, excluding the heat conduction mechanism through the 

(35) 

liquid, One can now write the heat balance equation as 

+ PBH - KgGg (36) 

where AT is the effective temperature difference between the melt and its 
surrounding environment. A comparison of equations (35 a,b) and (36) 
shows that the term h^AT must equal about 4 cal/cm^-sec and -6.79 cal/ cm -sec 
for the above cases of silver eind germanium respectively. This is consis- 
tent with the idea that convection carries heat through the melt to the 
interface for the silver case and away from the interface for the ge’*- 
manium case. Figure 32 indicates the conditions that may lead to convec- 
tion in the two systems. 

By measuring the temperature at the ampoule wall and at the center 
of the sample it was determined that the isotherms for silver were convex 
(viewed from the hot zone) in the melt immediately ahead of the interface 
and concave for germanium. For silver this means that the walls were 
hotter than the center leading to a lower density near the wall and a 
higher density at the center. Under the influence of gravity, this 
radial density gradient will set up a convective cell in which the mass 
at the center fails while the mass near the waJ I rises. In the case of 
germanium, however, the opposite condition is found; namely, near the 
melt-solid interface, the center is hotter than the wall and hence the 
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quartz capillary ^ 

and thermocouple Qp pQOR QUAUTY 

assembly 





Figure 31: A sketch showing where temperature measurements are 

made relative to the isotherms. 


Silver (convex) 


Germanium (concave) 



Heat flow for the 
curved isotherms 


Mass flow for the 
curved isotherms 


Figure 32: Possible heat and mass flow pattern for the indicated 

interface shapes. 
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mass at the wall falls while that at the center rises. The convective 
flow direction for germanium will therefore be opposite to that of silver. 
Note that equations (31) and (36) can be solved simultaneously to deter- 
mine the necessary expression for the correction factor to be applied to 
the thermal conductivity ratio determined in this work. According to 
the measured and literature values of the thermal conductivity ratios, 
of Table 2, similar percent errors were found when the conductivity ratio 
was expressed as a number greater than unity; that is, the conductivity 
value in the numerator is that of the phase with the higher conductivity. 
Accordingly, for silver and germanium respectively, this procedure will 
give 





(Ag) 


(37a) 


and 




(Ge) 


(37b) 


where Q » pRH + 

The terms Q/G^K^ and Q/Gj^K^ are the respective correction factors 
which empirically account for the percent errors listed in Table 2 for 
the thermal conductivity ratios. The result that the errors listed in 
Table 2 for silver and for germanium are about Che same. 
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i.e., 33Z, is not understood and may be fortuitous, but suggests that 
there may be self-cancelling effects which lead to roughly the same 
percent error. As mentioned earlier, for lack of a better correction 
procedure, the one of assuming a fixed percent error will be adopted. 

B. Position Error 

The raw daca obtained in the melt-solid temperature measurements 
consisted of temperature versus time. Since the translated length, Z, 
is the product of the translation rate, V, and time, t, an expression can 
be written as follows 

Z - Vt 

so that 

dZ - Vdt + tdV (38) 

An error of 2 seconds is known to exist in the time, recorded by the pro- 
grammable thermometer used in this investigation. This follows from the 
fact that the meter spends 2 seconds on any channel before it reads 
another one. There is also an error of ± .036 cm/hr in the measured 
translate rate. This is actually the standard deviation of the trans- 
lation rate values measured for the motor speed of 2.78 cm/hr. Finally, 
it is known that each of the dynamic experiments took about one hour from 
the start time to the time of the interface temperature measurement. Using 
these values, equation (38) becomes 



dZ ■ [ 2.8 cm/hr x hr] + [l hr x ± .036/hr] ■ 0.37cm 

jOUU 


The Initial position of the bottom of the ampoule is set by eye judgement 
of a pointer on the furnace graduated scale, before each of the dynamic 
runs. This may be the largest source of absolute position error and is 

I 

Conservatively assigned an error of 0.10 cm. These two errors add up to 
a total of ± .137 cm for. the interface position determined in this work. 
This means that for a translation rate, 2.78 cm/hr, a positional change 
in the ampoule of (2.78 ± .137) cm is read instead of 2.78 cm after a 
period of one hour. This gives an error of about 5Z in the position 
readings. Note that the measured temperature is affected by only that 
part of the position error that is time dependent. In the above case then 
only the .037 cm error affects the temperature readings from the "Tem- 
perature vs. distance" curves. 

It has been observed in connection with figure 23 that the melt- 
solid interface position of Pb Sn Te was essentially constant. How- 

X^X X 

ever, note that from figure 27 the freezing temperature is lowered as 
freezing progresses. This means a slight shift in the freezing isotherm 
toward the cold zone. The freezing temperature changed by about 11°C 
between the first and last interface position monitored. From table 3 
the gradient in the melt of the Bb Sn Te is about 80°C. For this type 

X“X X 

of temperature gradient, an 11°C change in freezing temperature should 
correspond to a positional shift in interface of about 0.13 cm. While 
this is consistent with figure 23 which shows a slight shift of about 0.1 


In Interface positions monitored at the slow translation rate of .046 
cm/hr, this amount of shift is nevertheless within the .137 cm calculated 
for the error in the experimental determina "on technique of the inter- 
face position and hence will be neglected. 

C. Concentration Profile Error 

In order to check the precision of the energy dispersive x-ray 
analysis, a ntut!.er of concentration determinations were made at one point 
on a sample using the same amount of count time, the same voltage and 
the same beam current. The standard deviation of the measurements is 
+ .006 mole fraction which is taken to be the precision of the measure- 
ments . 

During the alloy preparation, by weighing of the components, it is 
known that the initial concentration, C^, of the melt column of the 
Pb Sn Te is .2 mole fraction of SnTe. However, from the composition 
analysis discussed in Chapter IV, the best fit to the concentration pro- 
files of figure 27 gave an initial concentration value which differed 
from .2 mole fraction of SnTe by about 5%. This is interpreted to mean 
that the accuracy of the composition analysis is about 5%. 
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CHAPTER VI 
DISCUSSION 

A. Furnace Characterization 
1. Introduction 

It is generally known that a furnace environment affects the quality 

of simgle crystals grown in the furnace. Some of the furnace environment 

4 

properties which have been shown to affect crystal quality are the axial 

3 " ‘ ■ 

and radial temperature gradient and the insulation thickness or length 

of the adiabatic zone (region without a direct control heat source) . A 
tickler heater (booster heater) was used in this investigation to affect 
the axial temperature gradient. The tickler heater is a direct heat 
source and when placed at the end of the adiabatic zone, in effect re- 
duces the length of the adiabatic zone of the furnace. The extra heat 
generated by the tickler heater should shift the isotherms in the furnace. 
Again, it is known that the furnace temperature profile is dependent on 
the properties ’ of th*> load in the furnace. Chang et al. have shown 
that the isotherms of their modelled Bridgman type furnace do shift as 
the Biot number of the load is changed. The Biot number, hS/K, is de- 
fined as the ratio of heat loss from the sample surface to heat conduc- 
tion, where h is the load's heat transfer coefficient, S is the load 
radius and K is the thermal conductivity of the load material. 

A considerable amount of effort was put into characterizing the 
empty furnace used in this investigation in order to understand the 
environment in which tie solidification experiments were performed. 
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While temperature profile determinations are easier to make In an 

empty furnace than In a loaded one, both do obey the same fundamental 

heat flow principles. So a knowledge of the empty furnace profile of a 

crystal growth furnace can point to the trend In profile that may be 

» 

expected during growth. For Instance, an elimination or reduction of 
radial temperature gradients In an empty furnace usually corresponds to 
a similar sltiiatlon In a loaded furnace and represents a good crystal 
growth condition for obtaining a planar Interface. In effect, the result 
of an empty furnace characterization may be seen as a first approximation 
of the heat flow pattern to be expected in the loaded furnace, and 
hence the following discussion is concerned with characterizing the tem- 
perature profile under various furnace conditions. 

2. Effect of Plugging the Furnace Ends 

The result of leaving both the hot and cold zone ends unplugged In a 
vertical Bridgman type furnace is shown in figure 10. Notice the enormous 
difference between the temperatures of the near-wall of the quartz liner 
and the center-line of the furnace for any given axial position. This 
indicates the presence of large radial gradients. On plugging the fur- 
nace ends, an equilibration Is reached after about 3 minutes and de- 
pending on how effective the plugging is, the temperature differential 
between the quartz near-wall and the center line is either reduced or 
essentially eliminated, as shown in figures 11 and 15 respectively. 

Again, this Indicates a reduction or elimination of radial temperature 
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gradients. The temperature difference between the quartz near-wall and 
the furnace center-line may be caused by a "chimney" effect. This simply 
means that air is blowing into the furnace from the cold zone as hot air 
rises. However, due to frictional forces, the air near the wall is not 
rising as fast as the air in the center. So the wall retains its heat 
better than the center which is subjected to a flux of cooler air than 
that of the wall Also, the heat lost through the air flux will lower 
the temperature f both the hot and 'cold zones leading to a low axial 
gradient. If a tight plugging of the furnace ends is achieved such that 
essentially no air enters or leaves the furnace, then the radial gradient 
can be reduced remarkably and the axial temperature gradient increases. 

It has been mentioned that both the axia] nd radial temperature 
gradients affect the quality of crystals groim in the furnace. Since 
plugging of the furnace ends affects these furnace properties, as dis- 
cussed above, it must invariably affect the quality of crystals grown 
in the furnace. Therefore, proper plugging is essential to crystal 
growth in a vertical Bridgman type furnace in order to avoid the material 
problems due to high radial and low axial temperature gradients. Im- 
proper plugging leads to air flow in the furnace which creates a con- 
dition for the furnace wall to be substanially hotter than the center. 

3. Effect of Loading 

The experimental results show that a load Introduced into the fur- 
nace will have a temperature profile different from that of the empty 
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furnace. From Table 1, It can be seen that for the same furnace and 
tickler heater settings, the empty furnace profile was changed from a 
gradient of 127“G^cm at 900’c toU8‘’C/cm, 99‘’c/cm and Ss^C/cm by Pb^_^Sn^Te, 
germanium and silver respectively. This shows that as the load in the 
furnace changed from a low thermal conductivity material, Pb Sn Te, to 

X 

a higher thermal conductivity material, silver, the amount of change in 
the temperature gradient of 900°C Increased. 

The ability of a load to conduct heat means that if placed in a tem- 
perature gradient, the load will tend to go toward a thermal equilibrium 
by trying to eliminate the temperature gradient. So the load placed in 
the Bridgman type furnace will spread heat from the hot to cold zone in 
an attempt to remove the imposed temperature gradient. However, since 
there is a constant heat source set to maintain the gradient, the effect 
of the load is to reduce the temperature gradient. The higher the thermal 
conductivity, the more effective the load in spreading heat and reducing 
the temperature gradient. 

Figure 24 shows that the loaded furnace profile continuously 

shifted as the freezing of the silver progressed. This continuous shift 

in profile may be related to a change in the value of the overall Biot 

nvunber seen by the various zones of the furnace. As mentioned in the 

2 

introduction of this section, it has been shown theoretically that the 
Biot number, hS/K, is a parameter which affects the furnace temperature 
profile. An increase in this parameter is shown to correspond to a 
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shift In the Isotherms toward the cold end of the furnace plus a decrease 
in the spacing between the isotherms. The decrease in the spacing of the 
isotherms corresponds to a higher axial temperature gradient. Clyne^ 
has shown for aluminum that the shift in the thermal profile of the load 
is dependent on the proportion of the lengths of the specimen that is 
located in the hot and cold zones. 

For the same material, say silver, the heat transfer coefficient, h, 
is not expected to change substantially going from the melt into the solid 
phase, as compared to a factor of about two change in thermal conductivity 
values recorded for most semiconductors and metals. So changes in the 
value of h is neglected. For silver with = 2K^, this means that as 
the solid is formed, the Biot number, hS/K, is reduced in the cold zone. 
Also, the proportion of the entire starting column of material that is 
in the cold zone is increasing leading to a better coupling with the 
cold zone. These two conditions of a high solid thermal conductivity 
and a longer length of solid in the cold zone reinforces each other and 
leads to a shift in isotherms toward the hot zone as shown by figure Ik. 

The above discussions seem to indicate that the thermal conductivity 
of the load in a furnace is a factor that affects how the load tempera- 
ture profile compares with the empty furnace profile. The higher the 
thermal conductivity, the more the deviation from the empty furnace 
profile. As the load is traversed through the furnace, t' e change in 
the proportion of the length of the load that is in the cold zone causes 
a shift of isotherms toward the hot zone. 
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4» Effect of Tickler Heater 

In the furnace configuration used In this Investigation, the tickler 
heater Is placed within the adiabatic zone starting at the boundary of 
the hot and the adiabatic zones. As can be seen from figure 11, this Is 
within the region In which the furnace temperature profile rolls off from 
the hot zone to the cold zone. As already mentioned, the tickler heater 
In this furnace configuration basically reduces the length of the adia- 
batic zone. The heat generated In this region not only results In a 
positional delay In the roll-off "knee" but also Increases Its tempera- 
ture. This c.an be seen by comparing figures U and 15. The tickler 
heater generates hear only on the hot zone side of the insulation or 
adiabatic zone; whereas the temperature on the cold zone side Is fixed 
by the control thermocouple In the cold zone. Since as discussed above, 
the roll-off from the hot to the cold zone is delayed, it means that the 
gap within which the profile has to roll-off and fall to the fixed cold 
side temperature is smaller leading to a higher axial temperature gradient. 
As will be discussed in the next section, the above mentioned condition 
corresponds to a shift in isotherms toward the cold zone. 

B. Melt-Solid Interface Positions 

] . Introduction 

As discussed In the preceding section, a load put into a furnace 
can alter the temperature profile of the furnace. However, the dis- 
cussion did noL include the effect of the movement of the load, which 
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will be discussed here. Fu and Wilcox have shown theoretically that 

the furnace isotherns shift toward the cold zone as the Fcclet number , PCt 

increases. The Peclet number, VpC S/K, is defined as the ratio of heat 

P 

carried by motion of the load to heat conduction, where V is the load 
translation rate, p is the density, is the heat capacity, S is the load 
radius and K is the thental conductivity. 

The tickler heater has been shown to generate heat within the roll- 
off region, thereby increasing the axial temperature gradient of the 
empty furnace. In this section, this result will be correlated with 
shifts in the melt-solid interface and increases in the temperature 
gradient in the melt Immediately ahead of the interface. These are im- 
portant parameters In crystal growth by the Bridgman technique. Finally, 
the materials constants which were the offshots of the technique used 
in the Interface position determination will be discussed. 

2 . Tickler Heater Effect on Interface Positions 

The melt-solid interface has been previously referred to as repre- 
senting the isotherm that corresponds to the freezii'nj temperature of the 
material on the assumption that there is no compositional variation 
across the interface. In this sense then, any condition that shifts the 
furnace isotherms also shifts the interface. 

It has been shown In the preceding section that the extra heat gene- 
rated by the tickler heater raises the temperature of the roll-off "knee" and 
shifts the profile toward the cold zone. This shift means that any 
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given axial position within the roll-off region sees a higher temperature 
than It did without the tickler heater. This Is another way of saying 
that the position of a given Isotherm, 900°C for example. Is taken by a 
higher temperature Isotherm such as 950° while the 900°C Isotherm Is dis- 
placed further toward the cold zone. The Interface, as an Isotherm, also 
gets shifted for the same reasons. Of course, the longer the tickler 
heater length and/or the larger the tickler heat temperature setting, the 
more heat is generated In the 'oll-off region and the more the shift In 
Interface position. The heat generated by the *:ickler heater has also 
been shown In the preceding section to result in an Increase In the empty 
furnace axial temperature gradient. Recall that the empty f 'nace profile 
points to the trend In the profile to be expected during ;he crystal 
growth. So It Is expected .hat for the sane reasons given In the section 
on furnace characterization that the temperature gradient in the melt of 
a given material should increase with tickler heater length and/or setting. 
This is consistent with the results in Table 3 . 

3. Effect of Translation Rate 

2 

It has been pointed out in the introduction that a theoretical work 

by Fu and Wilcox has shown that as the Pec let number , Pe, increases, the 

furnace isotherms and therefore the melt-solid interface are shifted 

toward the cold zone. Note that from the definition of Pe, VpC S/K, for 

P 

a given material, the only variable is the translation rate, V. The 
above theoretical result then corresponds to a shift in Isotherms with 
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an Increase in the load translation rate. Also, equations (15 a,b) can 
be used to show that the dimensionless temperature, get bigger as Pe 
gets bigger. Again, for the jame material, this means that a given iso> 
therm Is replaced by a higher temperature one as the ampoule translation 
rate Is Increased. 

The above mentioned theoretical results are consistent with figure 

23 which shows that the melt-solid Interface of the same material 

Pb. Sn Te shifted toward the cold zone with an Increase In the transla- 
1-x X 

tlon rate, V, from .046 cm/m±n to .178 cm/min. Physically, the thermal 
mass of the load may be responsible for the load not d^'-fs gating heat 
as fast as it is being moved. So some of the heat carried from the hot- 
ter reglcr. : Is moved to the cooler region. This will make the cooler 
region hotter than it should normally be in a stationary condition. 

Of course, the higher the rate of sample translation, the less time the 
load has to dissipate heat to the immediate environment, the more heat is 
carried down to the cooler sections and therefore the further down the 
position of a given isotherm. 

4, Materials Constants 

The freezing point measured for germanium in this work is 937.3 ± 

1.5°C which is in agreement with the current literature value^^ of 937. 4°C. 
For silver, the measured . ’eezing point, 960.1 ± .9°C, is about 2°C 
lower than the current literature value of 961.93°C. These results may 
be related to the purity of the samples. The germanium used in the 
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Invaatigation was 99.9999 parcant pure wharaas the silver was only 99.99 
percent pure. From the phase diagrams of silver and oxygen or silver and 
sulphur, it can be seen that a small amount of these elements can sup- 
press the freezing point of silver significantly. In fact, an oxygen 
content of about 0.0001 mole fraction, according to the phase diagram, 
will suppress the freezing point of silver oy about l^^C. So the slightly 
low freezing point measured for silver is attributed to impurities like 
oxygen in the sample. 

Another material constant determined in this investigation is the 
melt-solid thermal conductivity ratios for the samples. The results for 
germanium and silver show an error of about 33X when compared with 
selected literature values, as shown in Table 3. This error has been 
discussed extensively in Chapter V. The total error is likely a combi- 
nation of errors due to the presence of the quartz capillary tube, the 
curvature of the interface and convection in the melt caused by radial 
temperature gradients. By taking temperature rsadings at both the wall 
and center of samples simultaneously, it was determined that the average 
radial gradient could be as high as 20*‘^C/cm which suggests that locally 
the radial gradient may be significantly higher. Hie radial temperature 
gradient creates a radial density gradient which coupled with gravity 
c.an Initiate and maintain convection in the melt. The calculations of 
Chapter V showed that the largest errar is likely that resulting from 


convection. 
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A generally accepted assumption in single crystal growth by the 
Brldgman-Stockbarger technique Is that solid-state dlffus^lon Is negli- 
gible compared to diffusion In the melt. As a result, mass redistribution 
during crystal growth from the melt Is explained by mass transport mecha- 
nisms In the liquid only. In the growth of Pb, Sn Te, the pseudobinary 
phase diagram of figure 8b suggests that a tln-rlch layer should build 
up In the melt In front of the Interface since the equilibrium segrega- 
tion coefficient,^^ k. Is less than unity. A number of things can happen 
to this tln-rlch layer. It can be completely or partially swept away and 
homogenized with the bulk melt by convection, Ir may remain undisturbed In 
the absence of convection. Of course, each of these conditions should 
yield a different final solute redistribution curve. 

One possibility Is the steady-state case In which there is no con- 
vection In the melt, and liquid mass transport Is by diffusion only. 

In this case, according to equation (7a), one expects a concentration 
profile that should have a constant composition over most of the length 
of the grown rod. Obviously, the data points in figure 27 do not exhibit 
any appreciable constant composition portion and so do not fit the steady 
state theory. Two other solute redistribtuion theories are; (a) 
complete mixing in melt by Pann^^ and (b) partial mixing in the melt 
with a diffusion boundary layer of thickness, 6, by Burton, Prim and 
Slichter^^. Both theories describe the solute redistribution according 


to equation (12) 
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except that in the latter theory, the equilibrium distribution coeffi- 
cient, k, is replaced by an effective distribution coefficient, 
equation (26) . Note that this partial mixing case is growth race depen- 
dent, through the definition of given by equation (10) 


k 


eff 


k 

k+(l - k)e 


-R6/Dl 


( 10 ) 


where R * crystal growth rate 

6 = effective diffusion layer thickness 

= diffusion coefficient in the melt 

37 

It is known from literature that the equilibrium segregation coef- 
ficient, k, ranges from .60 to .74 for liquid composition of about .2 to 
.5 mole fraction of Snie. However, when the data of figure 27 were 
fitted to equation (12), k values of 0.79 and 0.81 were obtained for the 
respective growth rates of .046 cm/min and .178 cm/min, both of which 
differed from the expected range values for the equilibrium segregation 
coefficient. This suggests that of the three theories discussed above, 
the BPS boundary layer approach gives the best description of the data 
measured in this investigation. 

From the above discussion, one is lead to believe that figure 27 
is indicative of the presence of free convection in the melt, with a 
boundary layer of thickness, 6 . While the difference between the 
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measured k and the literature value of k can be attributed to other 
err 

causes such as a breakdown In the assimptlons leading to equation (26), 

It is nevertheless reasonable to expect free convections to result from 

the radial temperature gradient which is known to exist in the melt and 

hence give rise to the k^^^ value of the BPS theory. 

From equation (10), a characteristic value, <S/0j^, for systems with 

free convection present in the melt can be calculated from a knowledge 

37 

of k , R and For the case of lead-tin telluride, the selected 

k value is .69. For the growth rate of .046 cm/min and .178 cm/min, 
the respective 6/D values are calculated to be 685 sec/ cm and 212 sec/cm. 

Lt 

A greater confidence is placed on the value of 685 sec/cm which is for 

the slower growth rate because it is based on many data points as opposed 

to that for the fast growth rate. This fact is evident in figure 27. 

38 

Clark has determined the diffusion coefficient, D , of SnTe in PbTe 

Lj 

as 7x10“^ cm^/sec. Using this D value, the thickness of the diffusion 

Ij 

boundary layer is calculated to be about .05 cm which is the same order 

35 

or magnitude as the 6 values listed in the literature for alloy systems 
undergoing progressive freezing with free convection present. 

C. Growth Rate 

The results shown in figure 25 indicate that the growth race of silver 
varied as freezing progressed and also that at all times it remained 
higher than the ampoule translation rate. In contrast, the growth rate 
of germanium and that of lead-tin-telluride remained approximately constant 
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and equal to the translation rate. As already mentioned in Section A of 
this chapter, the profile s^en by the load in a furnace seems to be de-> 
pendent on the proportion of the load lengths that are in the hot or cold 
zone. As more of the load enters the cold zone, there is a better 
coupling between the specimen and the cold zone while at the same tine 
there is a weak coupling with the hot zone. The net effect is that a 
given isotherm, e.g., the freezing isotherm, migrates tov;ard the hot zone. 
However, for a poor conductor like the semiconductors, germanium and lead- 
tin-telluride, once the first portion is frozen, heat transfer to the 
cold zone via conduction through the specimen is poor or essentially 
eliminated. Because of its poor conductivity, the solid essentially 
acts as an “insulator", separating the melt-solid interface from the 
cold zone. In this case, coupling with the cold zone is poor. Therefore, 
the loss of heat necessary for continued growth is sustained only by the 
mechanical drop of the sampJ.e toward the cold zone. Tliis nay explain 
the result chat the growth rate of the semiconductors was equal to the 
ampoule translation or drop rate. For Che metal silver, on tho other 
hand, the better coupling through the high thermal conductivity solid 
phase sustains a continued migration of the interface toward the liot 
zone, i.e., a continued displacement of the freezing isotherm Coward the 
hot zone. This added heat loss through the growth crystal is seen (from 
figure 25) to force the freezing isotherm to migrate faster than the 
ampoule drop rate and hence the continued change in interface position. 
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Figures 25 and 26 also show that the rate of migration of the Inter- 
face slowed down during the solidification process until about 66% of 
the melt column was frozen and then the migration of the Interface slowed 
a speed-up in the rate of advancement with translation of the boult . This 
may be related to the nonsymmetric nature of the furnace temperature pro- 
file, which goes from a low through a high and back to a low magnitude 
of the temperature gradient. Because the gradients were only measured 
near the Interface, there is no gradient data from which one could cal- 
culate the expected isotherm shift with translation rate and boule 
length. However, qualitatively it is possible for the boule to cause 
isotherm shifts as varying rates as a function of the translated boule 
distance, z. 

The resulLs suggest that for materials with very low thermal con- 
ductivity, the only means of effective heat loss that can ccntribute to 
solidification rate is provided by the sariple drop rate into the cold 
zone, However, for materials with high thermal conductivity, there is 
the .ihift in freezing isotherms due to better coupling with the cold 
zone. This plus the ampoule drop rate may make the rate of interface 
movement or the growth rate higher than the sample translation rate. 
Obviously, the magnitude of this effect is affected by the translation 
rate. 

D. Solute Redistribution 

The changing interfacial temperatures of the Pb, Sn Te due to 

1-x X 

changing compositions were determined from the melt-solid temperature 
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measurements. Examples of the measured freezing temperatures are 
shown In figure 27. From this figure, then both the concentration and 
the freezing temperature corresponding to that concentration can be 
determined at the chosen points. In effect, the solidus curve of the 
material can be deteoilned from the combination of x-ray analysis, com- 
position data, and the melt-solid temperature measurements technique 
used in this investigation. However, the pairs of concentration and 
temperature data taken from figure 27 tend to fall on the llquldus 
line of the pseudobinary phase diagram of figure 8b.. The samples used 
for the x-ray analysis were obtained by cutting across the rod crystal 
at the chosen points along the axis. This was then prepared as dis- 
cussed In Chapter III. The x-ray analysis was then done at points on 
the cross sectional surface on the assumption that the crystal was 
grown with a planar Interface. However, as already mentioned in this 
chapter, temperature measurement's made simultaneously at the wall of the 
ampoule and the center of the sample indicate the presence of radial 
temperature gradients which consequently correspond to a curved isotherm 
and hence a curved melt- solid interface. For the case of lead-tin- 
telluride, the isotherms and hence the melt-solid interface are concave 
(viewed from the hot zone) . The amount of curvature is represented by 
the magnitude of the radial temperature gradient. Unfortunately, the 
i;adial gradients were not measured for this case. In order to estimate 
the radial gradient, a proportionality between the axial temperature 
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gradient of lead-tln-tellurlde and its radial gradient will be assumed to 
be the same as that measured for germanium. The measured axial tempera- 
ture gradient for germanium is 99®C/cra from Table 1. A maximum radial 
gradient of 20°C/cm was also measured for germanium at the temperature, 
900°C. The axial temperature gradient for lead-tin-telluride is from 
Table 1, 118°C/cm. From this data, the radial temperature gradient of 
lead-tin-telluride is assumed to be 24°C/cm. 

Figure 30 has already been used to describe the radial positions at 
which the temperature measurements were made with reference to the curved- 
interface. Now if the sample is cuu in a planar cross section and x-ray 
analysis performed at the center, it becomes obvious that the measured 
concentration will correspond to a different isotherm than that at the 
point where the temperature reading was taken. The center of the x-ray 
samples is about 3.5 mm from the radial position where the temperature 
readings were taken. For a radial gradient of 24°C/cm, this center posi- 
tion therefore corresponds to an isotherm 8.4°C lower than the measured 
one. In order to obtain an estimate of the correct solidus temperatures 
at the radial points, the compositions were measured by x-ray analysis, 
one has to correct the measured temperature by subtraction of 8.4°C. Table 
7 shows the corrected solidus temperature which are to be compared with 
those values taken from the bulk phase diagram of figure 8b for the 
selected measured concentrations. The corrected solidus temperatures 
compare very well with those from the phase diagram. IJhile one can 



Measured 
Concentration 
Mole Fraction SnTe 

Measured Solidus 
Temperature Cor- 
rected for Interface 
Curvature 

Solidus Temperature 
Taken from Reference 
30 or Figure 8b 

.186 

901 °C 

098°C 

.191 

897°C 

396°C 

.206 

895“C 

894®C 

.220 

893°C 

892. 3°C 

.242 

890 °C 

039. 2°C 


Table 7: 


Measured solidus ter.iperatu. e corrected for the melt-solid 
interface curvature. 
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use the melt-solid temperature measurement technique and the x-ray analy 
ses to generate the solidus line of a given alloy system, care must be 
taken to properly correlate the temperature measurements and the compo- 
sition measurement so that they correspond to the same Isotherm and 


Interface. 
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(MAPTER VII 

SUMMARY Ain) CONCLUSION 


A. Growth Rates 

The instantaneous growth rate of directionally solidified silver was 
shown to remain higher than the tunopule translation rate all through the 
growth period monitored. The growth rate was continuously changing all 
through the growth period and at the early stages of growth was as much 
as 110/8: higher than the ampoule translation rate. On the other hand, 
the growth ri\tes of both Ge and Pb. Sn Te were essentially constant and 

X 

equal to the ampoule tmaslatlon rate within the monitored growth period. 

The continued migration toward the hot zone of the freezing isotherm 
of silver can be related to better heat coupling wi-.h the cold zone as 
more material is solidified. The apparent rate increase after about 66^ 
of the melt column has been frozen is possibly a result of the non- 
symmetric nature of the furnace temperature profile. 

From the discussion of the results, it is concluded that the growth 
rate of poor conductors like Ge and Pb Sn Te is determined by the sam- 
pie drop rate over the range of growth rates studied (.046 cm/min - 
.178 cm/min). However, the growth rate of high conductivity materials 
like metals may differ from the sample drop rate because of the added 
shift in the freezing isotherm due to better heat coupling with the 
cold zone . 

B. Materials Constants 

The measurement of temperature profiles on either side of a solid- 
liquid interface is an experimental technique employed in this research 
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to provide freezing temperatures of a 99.9999% pure germanium sample and 
a 99.99% pure sliver sample. The respective measured freezing points 
were 937.3 ± 1.5°C and 960.1 ± 0.9°C, which are In good agreement with 
literature values^^ of 937. 4°C and 961. 8®C. The discrepancy observed 
for silver Is consistent with a small Inpurlty level (e.g., .01% oxygen). 
This experimental technique was also used to determine the freezing 
temperatures of an alloy system at recorded points along the rod, as 
freezing progressed. Using a quantitative x-ray analysis, the composi- 
tion of the points for which freezing temperatures were recorded can be 
obtained, thereby completing a temperature versus composition curve which 
should correspond to the solidia line of the alloy system. A comparison 
of the composition-distance profile for lead-tin-telluride showed that of 
the three theoretical curves available In the literature, the one of Burton 
et al.^*^ (BPS theory) most closely fit the data. The BPS theory assumes 
a boundary zone of width 6 Into the melt through which mass transfer 
is by diffusion. Partial mixing of the melt and negligible diffusion 
in the solid is assumed. This theory leads to a continuously varying 
composition along the grown rod as shown in figure 27. Based on this 
theory, an effective segregation coefficient, of .790 and .810 

for the Pb, Sn Te was measured for the translation rates of .046 cm/min 
1-x X 

nad .178 cm/min respectively. 

38 

Using a diffusion coefficient calculated by Clark , the thickness, 

6, of the diffusion layer was calculated in this Investigation for 
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Pb, Sn Te to be 0.047 cm and 0.015 cm for the translation rates of 
1-x X 

0.046 cm/jidn and 0.178 cm/min respectively. The value of 0.047 cm is 
taken to be the more reliable value because it is based on a larger 
set of data points. 

Another direct measurement from the temperature profiles is the 

melt-solid thermal conductivity ratio of Pb. Sn Te. This is determined 

to be 2.33 ± .05. The solid thermal conductivity. Kg, of Pb^^^^Sn^Te at 

36 

the freezing point (x ■ .2) is given in the literature as 0.023 
Watts cm”^ K”^. With this Kg value and the ratio given aoove the 

melt thermal conductivity, K^^ can be calculated to be about 0.054 Watts 
cm“^ K Of course, this value is expected to vary with composition. 

C. Tickler Heater and PluRglng Effects 

A tickler heater (booster heater) has been used effectively to 
increase the temperature gradient, G , in the melt ahead of the inter- 

Li 

face. With the furnace configuration used in this investigation, the 
tickler heater essentially reduces the length of the adiabatic zone 
and generates heat within the temperature profile roll-off region of 
the furnace. In so doing, the temperature of the roll-off "knee' is 
elevated and the location of the entire roll-off region is displaced 
toward the cold zone of the furnace. This leads to an increase in axial 
temperature gradient, which has important implications for the growth of a 


single crystal. 
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The change In furnace temperature profile due to the tickler heater 
also corresponds to a condition In which any given Isotherm within the 
roll-off region Is pushed down toward the cold zone. This also Includes 
the freezing Isotherm or the melt-solid Interface which was shown to move 
toward che cold zone with an Increase In the extra heat generated In the 
roll-off region. The Increase In heat generation has been achieved by 
increasing the length or temperature setting or the tickler heater. 

Finally, it has been shown that improper plugging of the furnace 
ends may lead to a convective cold air flux in the furnace which makes 
the furnace liner walls hotter than the center line. This situation is 
shown to lead to high radial and low axial temperature gradients. There- 
fore, proper plugging of the furnace ends is essential to avoid to the 
materials problems caused by high radial and low axial temperature gra- 
dients . 

D. Load Material and Translation Rate Effects 

The temperature profile seen by the load during bulk cr>'stal growth 
using the Bridgman-Stockbarger technique is shown to change as more of 
the melt column is frozen. The data suggest that the temperature profile 
depends on the proportion of sample length in the hot and cold zones 
of the furnace and on the relative and absolute values of the thermal 
conductivities of the melt and solid. Figure 24 shows that for a high 
thermal conductivity material, silver, the profile sifts toward the hot 
zone. This shift in profile was not detectable for either Ge or Pb, Sn Te. 
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This result was interpreted to mean that if the thermal conductivity of 
the load is negligibly small, then the length of material in the hot or 
cold zone has a negligible effect on the furnace temperature profile. 

Table 1 shows how the furnace profile can be flattened to low 
axial gradient as the load material is changed from a poor conductor, 

Pb^ ^Sn^Te to a good conductor, Ag. The lowering of the axial tem- 
perature gradient corresponds to a shift in the isotherms toward the hot 
zone, and also an increase in the spacing between isotherms. However, 
the translation of the load through the furnace shows the opposite ef- 
fect of pushing isotherms toward the cold zone. This is demonstrated 
by the shift in freezing isotherms or the interface toward the cold 


zone with an increase in translation rate. 
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Data Reduction Program 
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A 


oRtGJuAi; r* - : ^ 

OF POOR QUALITY 


Listed belov/ is the program used in the reduction of the data from 
temperature versus time to temperature or temperature gradient versus 
distance. 


Ifc CREATE "TEMP4*',6 
2f CREAYE "GRADA",A 
30 CP ION BASE 0 
A0 ^SSIGh TO "TEMPA",! 

50 ASSIGN 08 TO "GRADA'^G 
60 INPUT "INITIAL DISTANCE*" , SI 
7u INPUT "INITIAL TEMP=",Y1 
8« INPUT "MOTOR SPEED =",C 
90 PRINT "MOTOR SPEED = ";C 
100 N»1 

110 PRINT "NO" , "TIME" , "DISTANCE" , "TEMP" , "TEMP GRAD" 

120 INPUT "TIME,TEMP = " ,X,Y 

130 S»S1+X*C 

140 IF SOS1 THEN 180 

1 50 PRINT 1 ,X ,S,Y ,0 

160 Z*S 

170 GOTO 220 

180 N»N+1 

190 D=(Y1 -Y)/(S1 -S) 

200 Z = S1 +X-ttC/2 

210 PRINT N,X,S,Y.Z,n 

220 S1=S 

230 Y1=Y 

240 PRINT «7;Y,S 

250 PRINT «8;D,2 

260 GOTO 120 

270 END 
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OF POOR QUALITY 


B 


The lattice constant calculation using values determined by x-ray 
diffraction (see Table 6A) is calculated by 

2 sine 

CuK - 1.5418 
a 

Using e - 27.67°, d - 3.224X 

For a cubic system (lead- tin- telluride has a Rocksalt structure) 

a^ « d^(h^ + + 1^) 

a - d*^ (h^ + k^ + 1^) 

From Table 6a, hkl « 200 


Therefore 

a - 6 . 42 X (Pb, Sn Te) 
1-x X 


From the ASTM cards 


^(PbTe) 

®(SnTe) 


6.443X 

6.302X 


Therefore 


^(PbSnTe) ’ ^(PbTe) " 


where x is the mole fraction of SnTe. 


X » - VbSnTe? . ,^^3 

.14 



appendix II 


Example of Melt-Solid Interface Position 
and Freezing Point Determination Curves 


TEMPERRTURE (°C) 
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OF POOR QUALITY 


MELT-SOLID TEMP. PROFILE (PbSnTe) 



0 .2 .4 ,B .0 I 1.2 

DISTRNCE (cm) 


Figure B1 . Second dynamic interface position and freezing point 
determination. Readings were tkaen every 30 seconds. 
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or!Gi:-;al i 

OF POOR QUALITY 

MELT-SOLID TEMP. PROFILE (PbSnTe) 


e .2 .4 .6 .8 1 1.2 



0 .2 .4 .6 .8 1 1.2 

DISTRNCE (cm) 


Figure 82. Third dynamic interface position and freezing point 
determination. Readings v/ere taken every 30 seconds. 
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OrilUlUAL » AA.- 

OF POOR QUALITY 

MELT-SOLID TEMP. PROFILE (PbSnTe) 


0 .2 .4 .B .0 1 1.2 



0 .2 .4 .6 .8 1 1.2 

□ISTHNCE (cm) 


Figure B3. Fourth dynamic interface position and freezing point 
determination. Readings were taken every 30 seconds. 
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OF POOR QUALITY 


MELT-SOLID TEMP. PROFILE (Ge) 



DISTANCE (cm) 

Figure B6. Second dynamic interface position and freezing point 
determination. 
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melt-solid temp. PP.OFILE (fig) 



Figure B8. First dynamic interface and freezing point determination. 
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